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Abstract 

1 

1/ 

A  steady  state  model  for  oxygen  atom  recombination  for  temperatures 
from  30C^K  to  200d?K  on  a  silicon  dioxide  surface  was  developed  based  on 
the  Langmulr-Rideal  heterogeneous  recombination  mechanism.  The  bonding 
of  atomic  oxygen  to  the  surface  and  the  thermal  desorption  of  atomic 
oxygen  from  the  surface  were  also  included  in  the  model.  The  hypothesis 
was  made  that  the  gas-phase  oxygen  atoms  combine  directly  with  the  oxygen 
atoms  that  constitute  the  silicon  dioxide  surface,  with  other  gas-phase 
oxygen  atoms  replacing  the  lost  atoms  that  were  bonded  on  the  silicon 
dioxide  surface  matrix.  The  model  agrees  with  the  experimental  data 
that  is  available  In  the  literature,  and  provides  an  insight  into  the 
processes  that  control  the  recombination  as  a  function  of  temperature. 

A  set  of  two-dimensional,  steady-state,  laminar  boundary-layer  calcula¬ 
tions  was  made  using  explicit  numerical  methods  to  demonstrate  the  use 
of  this  model  and  to  explore  the  rational  limits  which  may  be  placed  on 
the  role  of  catalysis  for  oxygen  recombination  on  a  silicon  dioxide 
surface  like  that  of  the  Space  Shuttle.  The  predicted  catalyclty  of  the 
silicon  dioxide  surface  did  reduce  the  overall  heating  below  that  which 
would  be  predicted  by  a  fully  catalytic  surface^ 


A  MODEL  FOR  OXYGEN  ATOM  RECOMBINATION 
ON  A  SILICON  DIOXIDE  SURFACE 

I.  Introduction 

Discrepancies  between  wind  tunnel  and  flight- test  heating  data  for 
the  reentry  phase  of  the  Space  Shuttle  flights  have  led  to  an  Interest 
in  catalytic  reactions  for  dissociated  air  species  on  silicon  dioxide 
(Ref  1-3).  In  particular,  detailed  klnetlc/physlochemical  models  for 
catalytic  recombination  of  nitrogen  and  oxygen  atoms  on  silicon  dioxide 
do  not  exist.  Models  for  these  species  are  important  because,  depending 
on  altitude/airspeed,  relatively  large  concentrations  of  atomic  oxygen 
and  nitrogen  are  expected  behind  the  strong  portion  (near-normal)  of  the 
bow  shock.  In  norequi librium  flow,  these  atomic  species  are  swept 
(convected)  along  the  heat-shield  surface  (Shuttle  tiles)  where  they 
kinetlcally  Interact  with  the  surface.  Since  the  recombination  process 
is  exothermic,  heating  rate  (and  resulting  surface  temperature)  predictions 
depend  heavily  on  the  extent  to  which  the  surface  is  catalytic.  High 
catalycity  leads  to  higher  heating  rates  and  higher  surface  temperatures 
and  low  catalycity  to  lower  heating  rates  and  lower  surface  temperatures 
(Ref  4-6).  Therefore,  surface  catalycity  provides  an  excellent  scapegoat 
to  blame  discrepancies  between  expected  heating  (developed  via  wind 
tunnel  testing)  and  experienced  heating  In  flight  testing. 

Although  such  discrepancies  may  well  be  due  to  wall  catalysis,  they 
may  also  be  due  to  other  nonequlllbrium  effects.  Scott  (Ref  7)  reviewed 
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the  nonequlllbrlum  calculation  techniques  developed  by  several  authors 
and  had  to  conclude  that  no  method  predicts  the  flight  measurements  of 
the  heating  over  the  entire  windward  centerline  surface  for  the  entire 
reentry.  He  noted  that  there  are  computational  difficulties  associated 
with  properly  characterizing  the  flow  field.  These  difficulties  arise 
from  the  complex  flow  and  the  gas-phase  chemical  reaction  kinetics;  in 
particular,  the  chemical  composition  of  the  flow  has  never  been  measured. 
There  are  disagreements  as  to  the  proper  treatment  of  the  boundary  condi¬ 
tions  at  the  surface  due  to  recombination  effects  (Ref  8).  Rosner  (Ref 
9)  points  out  that  the  recombination  data  that  has  been  inferred  from 
arc  jet  experiments  may  not  be  reliable  due  to  contaminants  In  the  arc 
jet  flow  field  (such  as  copper  from  the  electrodes)  and  since  flight 
conditions  are  not  accurately  simulated.  In  order  to  eliminate  discrep¬ 
ancies,  It  Is  Important  to  place  rational  limits  on  the  relative  role 
that  catalysis  can  play  in  the  overall  heating  problem.  Such  rational 
limits  cannot  come  from  comparisons  between  flight-test  and  wind  tunnel 
data,  as  apparently  has  been  done  for  the  Space  Shuttle  (Ref  10).  Such 
comparisons  simply  assign  a  catalyclty  which  makes  up  for  the  discrepancy. 
That  Is,  the  flight-test  data  are  due  to  thermal  energy  transfer  (Fourier 
conduction)  and  chemical  energy  transfer  due  to  recombination.  The 
perfect  gas  wind  tunnel  data  is  assumed  to  be  representative  of  equilibrium 
flow  with  a  fully  catalytic  wall  (Ref  7),  and  hence  a  catalycity  Is 
assigned  to  account  for  the  flight-test  heating  data  being  significantly 
less  than  the  wind  tunnel  data.  Rational  limits  for  catalysis  cannot  be 
made  In  this  manner.  However,  rational  limits  for  catalysis  can  be  made 
by  developing  models  which  derive  their  confidence  from  controlled 


experiments  specifically  tailored  to  provide  data  suitable  to  test  the 
validity  of  models  for  catalysis. 

Such  modeling,  verification  of  which  Is  based  on  controlled  experi¬ 
ments,  Is  the  object  of  the  present  research.  Because  of  the  extensive 
scope  of  the  most  general  task,  the  work  reported  here  deals  with  a 
subset  of  the  most  general  problem,  that  of  recombination  of  oxygen 
atoms  (only)  on  silicon  dioxide.  Application  of  the  results  of  this 
study,  then,  is  possible  for  a  limited  portion  of  the  Shuttle  reentry 
profile  when  oxygen  Is  dissociated  and  nitrogen  Is  not.  Because  of  the 
higher  dissociation  energy  for  nitrogen,  such  a  portion  of  the  reentry 
envelope  does  exist.  Thus,  a  partial  test  of  the  rational  limits  of  the 
role  of  catalytic  recombination  in  heating  predictions  for  the  Shuttle 
is  possible  since  the  coating  of  Shuttle  tiles  Is  predominantly  silicon 
dioxide  (Ref  11).  The  objective  of  this  research  Is  to  model  oxygen 
recombination  on  silicon  dioxide.  Second,  the  use  of  the  model  will  be 
demonstrated  for  the  purpose  of  showing  how  rational  limits  for  catalysis 
can  be  obtained. 
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II.  Modeling  of  Heterogeneous  Recombination 


Kinetic  studies  of  atom  recombl nation  on  most  surfaces  show  that 
the  reactions  are  first  order  processes  which  become  second  order 
processes  at  higher  temperatures  (c.f.,  below)  (Ref  12).  Yet  the  most 
common  method  of  collapsing  catalytic  data  uses  the  following  equation 


(Ref  13): 


Y  =  Yo  exP(-  UT> 


where  y0  and  E  are  constants.  yQ  Is  called  the  pre-exponential  factor, 
and  E  is  the  activation  energy.  R  Is  the  universal  gas  constant  that  Is 
used  when  energy  is  In  molar  units  and  often  is  replaced  with  the  Boltzmann 
constant  when  energy  Is  per  molecule.  Although  this  analytical  form  Is 
sometimes  useful  for  comparisons  with  other  recombination  coefficient 
data,  the  kinetic  details  are  lost.  The  "lumped  coefficient"  yq  tells 
little  about  what  it  Is  composed  of  in  terms  of  physical  and/or  chemical 
details.  Furthermore,  since,  as  will  be  shown,  this  form  requires  that 
the  reaction  be  first  order,  it  is  unable  to  describe  the  recombination 
coefficient  when  the  recombination  kinetics  are  second  order  or  when  the 
recombination  coefficient  goes  through  a  maximum. 

A  more  detailed  analytical  form  for  the  recombination  coefficient 
Is  given  by  Jumper,  et.  al.,  (Ref  14).  This  model  is  semi-classical  and 
describes  the  physical  and  chemical  microscopic  details  of  what  may  be 
occurring  on  the  surface  when  atoms  recombine  and  desorb  as  molecules. 

The  Jumper  model  allows  for  the  recombination  kinetics  to  transition  to 
second  order  at  elevated  temperatures  and  also  allows  the  recombination 


coefficient  to  pass  through  a  maximum  value  with  increasing  temperature. 
As  pointed  out  by  Jumper,  et.  al.,  in  Reference  14  for  fluorine  atom 
recombination  on  a  nickel  surface,  each  unknown  constant  In  the  recombina 
tlon  coefficient  model  led  to  a  specific  characteristic  of  the  data. 
Although  further  studies  are  required  to  verify  these  results.  Jumper's 
model  provides  much  more  Insight  Into  the  phenomena  possibly  occurring 
on  the  surface  where  atom  recombination  is  taking  place.  The  present 
oxygen  recombination  study  started  with  a  detailed  Investigation  of  this 
model.  In  the  following,  a  background  for  the  conventional  treatment  of 
catalysis  Is  addressed  followed  by  a  detailed  treatment  akin  to  the 
Jumper  method. 

Mechanisms  for  Heterogeneous  Recombination  of  Atoms 

Three  mechanisms  are  most  used  either  singularly  or  In  combination 
to  describe  atom  recombination  on  surfaces;  a  review  of  the  derivation 
of  the  so-called  Langmuir  adsorption  isotherm  Is  first  given  since  It  is 
useful  In  understanding  all  three  mechanisms.  The  Langmuir  adsorption 
isotherm  applies  to  the  Ideal  case  of  adsorption  where  there  are  no 
interactions  between  adsorbed  particles  (either  atoms  or  molecules). 

The  Isotherm  relates  to  adsorption  theory  as  the  perfect  gas  law  does  to 
the  study  of  real  gases.  The  adsorption  and  desorption  may  be 
represented  In  the  following  manner: 

k  A 

I  I 

A  +  -S-  — ^  -S- 


I 

where  -S-  represents  a  vacant  surface  site  and  A  Is  a  particle,  either 


an  atom  or  a  molecule.  If  we  let  0  be  the  fraction  of  the  surface  that 
is  covered  with  adsorbed  particles,  (1-0)  will  be  the  fraction  of  the 


surface  that  Is  bare.  The  rate  of  adsorption  will  then  be  k^p ( 1-©) » 
with  kj  being  a  rate  constant  and  p  the  pressure  of  the  particles  that 
are  being  adsorbed,  since  p  In  effect  represents  the  number  of  particle- 
surface  collisions  per  unit  area-time.  The  rate  of  desorption  is  then 
k_j0.  At  equilibrium,  the  adsorption  and  desorption  rates  are  equal 
leading  to  the  so-called  Langmuir  isotherm: 

0= - & —  (2) 

1  +  Kp 

where  K  Is  defined  as  k1/k_1  (Ref  12). 

Now  that  we  have  an  expression  for  the  fractional  surface  concentra¬ 
tion  0  in  terms  of  the  equilibrium  constant  K  and  the  pressure  of  the 
adsorbing  species,  we  can  look  at  the  three  mechanisms  which  are  normally 
used  to  model  the  heterogeneous  recombination  of  atoms. 


Langmulr-HInshelwood  Mechanism.  This  mechanism  Is  the  recombination 
of  two  atoms  adsorbed  on  neighboring  surface  sites  and  may  be  represented 
as: 


adsorption 

•  '  _ w 

A  +  A  +  -S-S-  ■ 

desorption 


^  A  recombination 

,  ,  desorption  „  ,  , 

-S-S-  -S-S-  +  a2 


the  recombination  desorption  rate  Is  proportional  to  0  times  0  since 
the  two  atoms  must  be  adsorbed  on  neighboring  sites.  Using  k  as  a 
proportionality  constant: 


Substituting  the  expression  for  0  from  the  Langmuir  adsorption  isotherm 
we  obtain: 


rlh  -  M — ^ — )2  (4) 

LM  1  +  Kp 

At  low  temperatures  the  surface  Is  most  likely  fully  covered,  0*1, 

which  means  that  Kp»l.  Therefore,  =  k  and  the  recombination  mechanism 

is  zeroeth  order  in  pressure,  l.e.,  It  does  not  depend  on  the  pressure. 

At  higher  temperatures,  the  surface  Is  sparsely  covered,  Kp«l,  and  R^ 

2  2 

=  kK  p  .  The  rate  goes  as  pressure  squared  and  thus  the  kinetics  are 
second  order  at  the  elevated  temperatures. 

Slow  Adsorption  Fol 1  owed  by  Rapid  Migration  and  Recombination  Mechanism. 
A  second  mechanism  Is  considered  by  Laldler  (Ref  15),  and  this  is  the 
surface  migration  of  chemisorbed  atoms.  This  mechanism  could  only  be 
significant  when  the  fractional  surface  concentration  0  Is  small,  i.e., 
the  surface  Is  sparsely  covered.  The  recombination  rate  would  be  propor¬ 
tional  to  the  adsorption  rate  and  given  as: 

Migration  =  kp(1  " 

If  0  1s  small,  the  recombination  rate  would  be  first  order  with  respect 
to  the  pressure  of  the  recombining  atomic  species.  However,  this 
mechanism  would  remain  first  order  even  at  higher  temperatures. 


Lanqmuir-Rideal  Mechanism.  The  third  mechanism  is  the  recombination 
of  an  atom  from  the  gas  phase  with  an  atom  that  is  adsorbed  to  a  surface 


site: 


adsorption  A 

<  •  _ »  > 

A  +  A  +  -S-S-  -  A  +  -S-S- 

desorption 


recombination 
desorption  „ 


-S-S-  +  a2 


The  recombination  desorption  rate  RLR  Is  now  proportional  to  0  times  p, 
the  pressure  of  the  atoms: 


Rlr  =  kp©  (6) 

Again  substituting  for  0: 

rlr  ■  ‘pfrH")  <7> 

At  low  temperatures  Kp»l,  and  the  kinetics  are  first  order.  At  elevated 
temperatures  where  Kp«l,  the  recombination  kinetics  are  second  order 
with  respect  to  the  pressure  of  the  atomic  species. 


Adopted  Mechanism 

In  order  to  allow  the  recombination  kinetics  to  transition  from 
first  to  second  order  with  respect  to  the  atomic  pressure  as  the  tempera 
ture  Increases  and  also  to  keep  the  model  from  becoming  unnecessarily 
complex,  the  Langmulr-Rldeal  recombination  mechanism  was  adopted  for  the 
model  for  oxygen  recombination  and  hydrogen  recombination  (see  Appendix 
A)  on  silicon  dioxide.  Although  this  mechanism  Is  acceptable  and  often 
used  to  explain  the  recombination  process  (Ref  12,  16-18),  it  is  not  the 


only  second  order  recombination  mechanism  at  elevated  temperatures. 

Wood  and  Wise  (Ref  19)  modeled  the  recombination  of  hydrogen  on  silicon 
dioxide  using  a  Langmulr-Rldeal  mechanism  from  300°K  to  500°K  and  a 
Langmulr-Hinshelwood  mechanism  above  500°K.  Hardy  and  Linnett  (Ref  18) 
point  out  that  both  the  Langmuir-Rideal  and  the  Langmulr-Hinshelwood 
mechanisms  require  the  thermal  desorption  of  atoms  to  be  large  at  the 
elevated  temperatures  where  the  kinetics  are  second  order,  and  therefore, 
that  the  order  of  the  recombination  process  will  not  exclude  one  of 
these  two  mechanisms  from  the  other. 


Components  of  the  Rate  Equation 

Let  us  now  look  at  the  microscopic  processes  which  are  possibly 
occurring  at  the  surface  upon  which  heterogeneous  reactions  are  taking 
place. 

Physical  and  Chemical  Adsorption.  At  a  surface,  a  molecule  or  atom 
may  experience  a  weak  interaction  due  to  Van  der  Waals  forces.  This 
weak  interaction  Is  physical  adsorption  and  Is  associated  with  llquefacti 
and  with  condensation.  If  a  molecule  or  atom  experiences  a  strong  Inter¬ 
action  at  the  surface,  the  particle  will  undergo  chemical  adsorption 
(chemisorption)  and  a  strong  bond  will  be  formed  with  the  surface.  In 
general,  the  heat  of  adsorption  for  physical  adsorption  Is  In  the  range 
2-15  kcal/mole,  whereas  for  chemisorption,  the  heat  of  adsorption  Is  In 
the  range  15-200  kcal/mole  (Ref  20).  Physical  adsorption  tends  to  occur 
only  at  temperatures  near  or  below  the  boiling  point  of  the  adsorbed 
particles  at  the  pressure  considered  (Ref  21).  Jumper,  et.  al.,  (Ref 
14),  show  that  physical  adsorption  Is  negligible  at  temperatures  above 


100°K  for  their  model  for  fluorine  recombination  on  a  nickel  surface. 


Their  soft  cube  model  provides  an  explanation  of  the  gas  particle-surface 
interaction.  For  oxygen  atom  recombination  on  oxide  surfaces  at 
temperatures  greater  than  30G°K,  the  only  significant  form  of  adsorption 
Is  chemisorption  (Ref  13).  Since  the  recombination  of  oxygen  atoms  on 
silicon  dioxide  from  300°K  to  2000°K  is  modeled,  only  chemisorption  was 
considered. 


Chemisorption  Rate.  The  rate  of  chemisorption  of  particles  (atoms 
or  molecules)  to  silicon  dioxide  can  be  expressed  as  follows: 

chemisorption  rate  =  NS  (8) 

where  N  is  the  surface  impingement  rate  per  unit  area  and  S  is  the 
sticking  coefficient.  The  surface  impingement  rate  is  assumed  to  be  gas 
kinetic  and  is  given  by  (Ref  22): 


N 


_  c 
n  T 


(9) 


where  n  is  the  number  density  of  particles  In  the  gas  and  c"  is  the 
average  velocity  of  the  particles.  The  average  velocity  is  derived  in 
kinetic  theory  as: 


8kBT  4 

c  *  ( — - — y 
'  -rrm  ' 


(10) 


where  kg  is  the  Boltzmann  constant,  T  is  the  absolute  temperature,  and  m 
is  the  mass  of  the  particles.  The  sticking  coefficient  may  be  expressed 
as: 

(ID 


S  =  so(i  -  0) 
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where  SQ  is  the  Initial  or  clean  surface  sticking  coefficient  and  9  is 
the  fractional  surface  concentration.  This  functional  dependence  upon 
SQ  and  0  assumes  (1)  that  any  particle  striking  an  occupied  site  is 
reflected  back  into  the  gas  phase,  (2)  each  site  can  be  occupied  by  only 
a  single  particle,  and  (3)  surface  mobility  of  chemisorbed  particles  is 
negligible. 

Thermal  Desorption  Rate.  The  rate  of  thermal  desorption  of  the 
particles  is  given  by  (Ref  23) 

thermal  desorption  rate  =59  (12) 

where  6  is  the  thermal  desorption  rate  per  unit  area  of  covered  surface. 
According  to  the  theory  of  absolute  reaction  rates  and  assuming  the 
ratio  of  the  activated  complex  partition  function  to  the  adsorbed 
particle  partition  function  is  unity: 


5  =  C  kBT  exp(- 
a  — c — 


where  C  Is  the  number  of  adsorption  surface  sites  per  unit  area,  h  Is 
a 

the  Planck  constant,  and  D  Is  the  thermal  desorption  energy. 

Recombination  Desorption  Rate.  The  rate  at  which  the  chemisorbed 
atoms  recombine  with  atoms  from  the  gas  phase  and  desorb  from  the 
surface  as  molecules  is  given  by 


Recombination  =  Surface 
Desorption  Impingement 

Rate  rate 
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which  is  analogous  to  the  collision  theory  results  for  homogeneous 
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reactions  (Ref  22).  The  four  factors  In  equation  (14)  may  be  written  as: 


•  C 

Recombination  =  N0  [exp  (-  t=t)]P  (15) 

Desorption  Rate  KB 

where  E  is  the  apparent  activation  energy  and  P  is  the  steric  factor 
which  must  be  empirically  determined. 

Surface  Catalyzed  Formation  Desorption  Rate.  If  a  molecule  is 
formed  as  a  result  of  a  gas  phase  particle  (atom  or  molecule)  striking  a 
chemisorbed  particle  (atom  or  molecule),  the  Langmuir-Rideal  mechanism 
can  again  be  used  to  describe  the  formation  of  product  molecules  and 
desorption: 


Surface  Catalyzed  .  F 

Formation  =  N0  [exp(-  -r=y)]P 

Desorption  Rate  KB 


(16) 


Rate  Eauatlon 


As  mentioned  In  the  previous  chapter,  air  consists  predominantly  of 
five  species  at  temperatures  below  8000° K  (Ref  22):  N,  N2,  0,  02,  and 
NO.  The  time  rate  of  change  of  the  number  of  chemisorbed  nitrogen  atoms 
per  unit  area  assuming  atoms  are  chemisorbing,  thermally  desorbing, 
recombining  via  the  Langmuir-Rideal  mechanism  and  subsequently  desorbing 
as  molecules,  and  forming  NO  via  a  Langmuir-Rideal  mechanism  and  then 
desorbing  can  now  be  expressed  as: 


dn 


Ut 


SN  =  S_  Nn  (1-9|0  -  SN0N  -  PNN  Nn  ^  exp(-  JJy) 
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(17) 
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where  n  is  the  number  of  chemisorbed  nitrogen  atoms  per  unit  area.  The 
SN 

subscript  N  on  the  parameters  signifies  that  the  variable  is  for  atomic 
nitrogen.  The  double  subscript  NN  and  ON  mean  that  the  variable  is  for 
the  nitrogen  recombination  and  NO,  nitric  oxide,  formation,  respectively. 
A  similar  equation  may  be  written  for  the  time  rate  of  change  of  the 
number  of  chemisorbed  oxygen  atoms  per  unit  area  when  oxygen  atoms  are 
undergoing  the  same  processes: 

dn  .  .  Enn 

-g  ■  So0  "o  O-V  -5o  0o  -  poo  No  Bo  T^T> 

■  PNO  "ll  ®0  exp("  1c^T!  (18) 


The  double  subscript  NO  signifies  the  surface  catalyzed  NO  formation  via 
the  collision  of  a  gas  phase  nitrogen  atom  with  a  chemisorbed  oxygen 
atom,  while  the  ON  subscript  in  equation  (17)  Is  for  the  formation  of 
nitric  oxide  via  the  collision  of  a  gas  phase  oxygen  atom  with  a  chemi¬ 
sorbed  nitrogen  atom. 

If  nitrogen,  oxygen,  and  nitric  oxide  molecules  are  assumed  to  be 
able  to  chemisorb  to  and  thermally  desorb  from  the  surface,  then  rate 
equations  may  be  written  for  the  time  rate  of  change  of  the  number  of 
chemisorbed  molecules  of  each  of  these  species  in  the  following  way: 


l=s0„  "n  »-eT>-Vl 


2  ■  s0„  "o'1  -  ®T>  -  so  ®0 
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Note  that  these  equations  (19)  -  (21)  have  only  one  means  whereby  the 

molecules  leave  the  surface  -  via  thermal  desorption.  0y  is  the 

total  fractional  surface  concentration  and  is  bounded  by  0  £  0y  <  1. 

0T  =  0N  +  0N2  +  G0  +  %  +  0NO  ^ 

This  expression  for  0y  may  be  substituted  into  equations  (17)  -  (21). 
In  steady  state,  these  equations  are  each  equal  to  zero,  and  therefore 
may  be  solved  for  the  fraction  of  the  surface  that  is  covered  by  each 
species: 


1  ■  °N2  ■  ®o  "  0O2  ■  9no 


1  +  _A_  +  PNN  exp(-  j-y)  +  P0N  N0  exp(-  ^l) 


So  nn 
°N  N 


1  “  °N  '  90  '  902  "  9N0 


s  Nn 
°N2  N2 


S  Nn 
°N  N 


1  "  ®N  "  GN2  '  902  "  9N0 


1  +  — +  P00  exp("  +  PN0  exp(‘ 


°0  0 


1  '  °N  ‘  GN2  "  90  "  0NO 


S  Nn 
°02  °2 


S  Nn 
°0  0 


mmmmm 


,  ■v  .  "s.  .  "V  \ 

"  % 


These  isolated  fractional  surface  concentrations  correspond  to  the  actual 
surface  coverages  when  the  other  species  are  not  present.  By  substituting 
these  definitions  for  the  isolated  fractional  surface  concentrations, 
equations  (28)  -  (32),  into  equations  (23)  -  (27),  we  obtain  the  following 
system  of  five  equations  in  terms  of  the  five  fractional  surface  concentra 
tions  for  the  five  species: 

1 

1 

W~o 
1 

1 

(33) 

This  system  of  equations  may  now  be  solved  in  a  straight  forward  manner 
for  the  fractional  surface  concentrations  using  Cramer's  Rule. 

Rather  than  give  the  details  of  this  solution.  It  is  sufficient 
here  to  comment  on  the  complexity  of  the  resulting  expressions.  The 
resulting  expressions  for  the  fractional  surface  concentrations  are  In 
terms  of  five  initial  sticking  coefficients,  five  thermal  desorption 
energies,  four  activation  energies,  and  four  steric  factors.  Some  of 


these  parameters  are  not  constants,  and  we  have  a  very  complex  system  to 
deal  with.  Considerable  further  simplifying  must  be  accomplished  before 
this  system  of  equations  can  be  handled  since  many  of  these  parameters 
are  unknown  and  must  be  experimentally  determined.  The  next  chapter 
discusses  the  simplifications  that  were  made  to  reduce  the  complexity  of 
the  solution  of  equation  (33). 


In  order  to  Unit  the  number  of  physical  processes  taking  place  on 
the  catalytic  silicon  dioxide  surface.  It  would  be  nice  to  consider  only 
the  0  to  Og  reaction.  Fortunately,  this  case  has  some  practical  signifi¬ 
cance.  During  the  latter  portion  of  the  reentry,  a  vehicle  like  the 
Space  Shuttle  has  decelerated  sufficiently  so  that  the  nitrogen  atom 
concentration  in  the  flow  over  the  vehicle  surface  Is  negligible  (Ref 
24).  Under  these  conditions,  air  consists  primarily  of  only  four  species 
N2>  0,  02,  and  NO.  Therefore,  the  equations  derived  In  Chapter  II  may 
be  simplified  to  a  case  which  can  both  be  handled  and  have  utility.  If 
we  further  assume  that  the  nitrogen  and  nitric  oxide  molecules  act  only 
as  diluents  and  as  such  do  not  chemisorb  to  the  surface  sites  on  the 
silicon  dioxide  surface,  we  are  left  with  only  oxygen  atoms  and  molecules 
actively  participating  on  the  surface. 

In  that  which  follows,  a  development  is  given  which  adapts  wholesale 
the  approach  of  Jumper,  et.  al.,  (Ref  14),  for  recombination  of  fluorine 
on  nickel.  As  will  be  discussed  in  Chapter  V,  this  approach  leads  to 
results  which  are  unphysical;  It  does,  however,  provide  a  starting  point. 
Using  the  approach,  then,  of  Jumper,  et.  al.,  (Ref  14),  and  the  above 
simplifications,  equation  (18)  reduces  to  the  following: 


We  can  solve  for  the  fractional  surface  concentrations: 
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Recombination  Coefficient 

The  recombination  coefficient,  y.  Is  defined  as  the  fraction  of 
collisions  with  the  surface  leading  to  recombination  (Ref  25).  Therefore, 
the  recombination  coefficient  Is  equal  to  two  times  the  recombination 
desorption  rate  divided  by  the  atom  impingement  rate.  The  factor  of  two 
accounts  for  the  fact  that  two  atoms  must  collide  with  the  surface  for 
each  recombination  with  the  first  atom  chemisorbing  to  the  surface.  The 
recombination  coefficient  Is  then  given  by: 

Y  =  2  P  0O  exp(-  yj)  (42) 

B 

If  we  substitute  equation  (40)  for  0, 


2  P  (9q)  [1  -  (0q  )  ]  exp(-  -p-y) 

u  o  u2  o  *B 


Alternatively,  In  terms  of  the  basic  parameters,  equation  (43)  Is: 
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In  either  form,  we  see  that  a  number  of  physical  parameters  are  needed 
to  proceed  further.  These  parameters  are  the  initial  sticking  coeffi¬ 


cients  for  atomic  oxygen  and  molecular  oxygen,  S_  and  ,  respectively; 

°0  °o2 

the  thermal  desorption  energies  for  atomic  oxygen  and  molecular  oxygen, 

Dq  and  Dq  ,  respectively;  and  also  the  sterlc  factor  and  the  activation 
energy  for  the  recombination  process. 

Equation  (44)  is  for  the  oxygen  atom  recombination  coefficient  when 
both  the  oxygen  atoms  and  the  oxygen  molecules  are  assumed  to  chemisorb 
to  the  silicon  dioxide  surface.  It  Is  Important  to  again  note  that  we 
are  assuming  that  both  oxygen  atoms  and  molecules  chemisorb  to  the  silicon 
dioxide  surface,  and  that  the  sticking  coefficient  is  given  be  equation 
(12),  where  SQ  may  be  a  function  of  temperature,  consistent  with  the 
model  for  fluorine  recombination  on  nickel  suggested  by  Jumper,  et.  al. 
Although  I  will  attempt  In  Chapter  V  to  argue  that  this  Is  not  the  case, 
it  Is  of  interest  to  explore  the  ramifications  of  these  assumptions  in 
light  of  Interpretations  of  the  effect  on  the  character  of  the  data  each 
influencing  parameter  has. 

Jumper,  et.  al.,  (Ref  14),  made  the  observation  that,  while  there 

was  an  Interrelationship  between  the  parameters,  each  had  a  distinct 

effect  on  the  character  of  the  data.  Equation  (43)  Is  similar  to  that 

of  Jumper's  model  for  fluorine  atom  recombination  when  the  temperature 

Is  sufficiently  high  so  that  physical  adsorption  Is  negligible,  with  the 

exception  that  Jumper,  et.  al.,  assumed  that  the  activation  energy  is 

negligible  for  fluorine  atom  recombination  on  nickel  and  thus  their 

equation-like  equation  (43)  does  not  contain  the  exp  (-  factor.  In 

B 

this  form,  it  is  clear  that  a  competition  for  surface  sites  between  the 


atoms  and  molecules  Is  indicated,  and  this  competitive  role  formed  the 
key  mechanism  to  Interpret  the  fluorine  recombination  data  by  providing 
the  ability  for  the  recombination  coefficient  to  first  Increase  with 
temperature  but  then  decrease  as  the  temperature  continues  to  Increase. 

According  to  the  explanation  of  Jumper,  et.  a!.,  at  temperatures 
where  only  chemisorption  is  significant,  the  recombination  coefficient 
increases  with  increasing  temperature  to  a  maximum  and  then  decreases  as 
the  temperature  Is  further  Increased.  The  recombination  coefficient  Is 
able  to  Increase  as  the  temperature  Increases  because  the  chemisorbed 
molecules  are  thermally  desorbing,  thereby  allowing  more  atoms  to  chemisorb 
to  the  surface  and  subsequently  recombine.  At  still  higher  temperatures, 
however,  the  chemisorbed  atoms  begin  to  thermally  desorb,  thus  allowing 
for  the  recombination  coefficient  to  go  through  a  maximum,  i.e.,  the 
thermal  desorption  Increases  with  Increasing  temperature  causing  a  decrease 
In  the  fractional  surface  concentration  of  the  atomic  species  and  a 
decrease  In  the  recombination  rate. 

Jumper,  et.  al.,  (Ref  14)  modeled  the  initial  sticking  coefficients 
for  both  atomic  and  molecular  fluorine  on  nickel  and  as  a  result  only 
the  sterlc  factor  and  the  two  thermal  desorption  energies  were  unknown. 

They  found  that  the  remaining  three  variables  affected  the  recombination 
data  In  distinctly  different  ways,  and  consequently  could  be  determined 
from  the  data.  There  Is,  of  course,  no  assurance  that  the  same  model 
that  worked  for  the  recombination  of  fluorine  on  nickel  should  work  for 
the  recombination  of  oxygen  on  silicon  dioxide,  but  both  the  model  and 
the  method  for  obtaining  the  missing  parameters  at  least  formed  a  starting 
point.  Indeed,  this  is  what  has  been  done.  Rather  than  describing  that 


IV.  Data  for  Oxygen  Recombination  on  Silicon  Dioxide 

Recombination  Data 

In  principle,  the  heterogeneous  recombination  of  atomic  oxygen  is 
an  elementary  surface  reaction.  However,  the  recombination  coefficient 
results  for  oxygen  recombination  on  silicon  dioxide  often  vary  by  a 
factor  of  five  within  a  single  laboratory,  and  by  up  to  two  orders  of 
magnitude  between  different  laboratories  (Ref  25).  Llnnett's  work  (Ref 
17,  26-28)  provides  the  recombination  coefficient  for  oxygen  recombination 
on  silicon  dioxide  by  two  different  methods.  These  results  are  considered 
the  best  data  to  date  and  have  even  been  used  for  calibration  purposes 
(Ref  9).  Figure  1  shows  the  recombination  coefficient  data  points  as  a 
function  of  the  reciprocal  of  the  absolute  temperature  as  taken  from 
Reference  17.  The  data  points  shown  are  for  one  series  of  experiments, 
but  the  flattening  of  the  curve  drawn  through  the  data  by  the  experi¬ 
mentalists  at  the  elevated  temperatures  was  confirmed  In  further  tests. 
Rosner's  work  as  reported  by  Breen,  et.  al.,  In  Reference  9  also  supports 
the  results  of  Linnett.  The  flattening  of  the  curve  In  Figure  1  appears 
to  be  due  to  a  maximum  value  for  the  recombination  coefficient.  It  is 
known  that  as  the  temperature  Is  further  Increased,  the  recombination 
coefficient  decreases  (Ref  9). 

As  shown  In  Figure  1,  Llnnett's  data  increases  in  value  nonlinearly 
as  the  temperature  Increases.  Greaves  and  Linnett  (Ref  17)  proposed 
that  the  activation  energy  Increases  from  1  kcal/mole  at  room  temperature 
to  13  kcal/mole  at  300°C. 


The  order  of  the  recombination  process  is  also  an  extremely  important 


piece  of  data.  The  recombination  process  is  first  order  with  respect  to 
the  oxygen  atom  pressure  from  room  temperature  to  350°K  (Ref  26).  The 
recombination  process  has  been  experimentally  verified  to  remain  first 
order  from  room  temperature  to  above  1100°K  (Ref  9).  As  was  described 
In  Chapter  II,  the  recombination  kinetics  on  most  surfaces  are  In  general 
first  order  processes  which  transition  to  second  order  at  higher  tempera¬ 
tures  (Ref  12).  While  no  empirical  evidence  exists  for  the  recombination 
mechanism  transitioning  to  second  order  for  oxygen  recombination  on 
silicon  dioxide.  In  view  of  the  above,  a  satisfactory  model  should  allow 
for  such  a  transition. 

Silicon  Dioxide  Surface  Structure 

In  each  of  the  various  forms  of  silicon  dioxide  the  silicon  atom  is 
bonded  in  a  covalent  manner  (Ref  29)  at  the  center  of  a  regular  tetra¬ 
hedron  (Ref  30).  The  tetrahedra  link  together  corner  to  corner  and  the 
different  structures  formed  determine  the  different  forms  of  silicon 
dioxide — cristobal ite,  quartz,  or  tridymite.  Figure  2  is  a  reproduction 
of  Wells'  Figure  288  and  shows  the  different  silicon  dioxide  structures. 

The  surface  structure  of  silicon  dioxide  must  satisfy  the  stoichio¬ 
metric  relationship  SIO2  In  addition  to  the  substructure  matrix  such 
that  each  silicon  atom  Is  bonded  to  four  oxygen  atoms  via  half  bonds: 
Si(i0)4.  That  Is,  each  oxygen  atom  in  the  substructure  is  shared  by  two 
silicon  atoms.  Heterogeneous  recombination  is  very  dependent  upon  the 
surface  structure  and  this  may  be  the  reason  for  such  large  differences 
in  experimental  results  (Ref  31).  Greaves  and  Linnett  (Ref  17)  discuss 


the  possible  surface  structure  when  hydrogen  atoms  are  present  In  the 
form  of  water  with  the  result  that  hydroxyl  groups  are  attached  to  the 
silicon  dioxide  surface.  They  mention  that  when  no  OH  groups  are  present, 
the  surface  structure  might  look  like  the  following: 


/  \ 
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where  the  oxygen  atom  that  has  the  double  bond  forms  the  surface  layer 
with  other  oxygen  atoms  that  have  double  bonds  to  silicon  atoms  below 
them. 


Removal  of  Surface  Oxygen  Atoms 

Experimental  evidence  that  is  independent  of  oxygen  recombination 
data  shows  that  oxygen  atoms  can  be  removed  from  a  silicon  dioxide 
surface.  Hickmott  (Ref  32)  studied  the  interaction  of  atomic  hydrogen 
with  glass  and  found  that  the  glass  provided  the  oxygen  atoms  which  were 
In  the  product  molecules. 

Campbell  and  Thrush  (Ref  33)  observed  BN  gas  downstream  of  their 
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Pyrex  flow  tube  following  nitrogen  recombination  studies.  Rahman  and 
Llnnett  (Ref  34)  state  that  nitrogen  atoms  evidently  bond  to  boron,  and 
therefore,  the  active  sites  may  be  the  boron  and  silicon  atoms  In  the 
glass. 

Evenson  and  Burch  (Ref  35)  found  that  the  recombination  coefficient 
for  nitrogen  atom  recombination  on  quartz  decreased  as  the  oxygen 
impurity  of  their  nitrogen  source  gas  increased.  That  Is,  the  oxygen 
atoms  act  as  a  poison  for  nitrogen  recombination  on  these  quartz 
surfaces. 

Rosner  (Ref  9)  used  X-Ray  Photoelectron  Spectroscopy  (XPS)  to 
analyze  silica  and  silicon  carbide  samples  exposed  to  an  oxygen  atom 
plus  nitrogen  atom  mixture.  Silica  samples  exposed  to  this  mixture  at 
1490°K  showed  that  nitrogen  atoms  reacted  with  the  silica  forming  either 
silicon  nitride  or  silicon  oxynitride  or  possibly  just  were  absorbed. 

When  silicon  carbide  samples  were  exposed  to  the  same  mixture  of  oxygen 
and  nitrogen  atoms  at  temperatures  up  to  1620°K,  a  surface  of  only  silica 
was  produced.  Oxidation  is  evidently  favored  over  nitrldatlon  at  tempera¬ 
tures  up  to  1620°K.  At  higher  temperatures  following  quenching  in  argon, 
the  silicon  carbide  samples  consisted  of  silicon  carbide,  silicon  nitrides, 
or  silicon  oxynitrides.  (This  surface  activity  will  be  referred  back  to 
in  Chapter  V  where  the  details  of  the  modified  model  are  given.)  An 
atom  that  is  chemisorbed  to  a  surface  is  essentially  the  same  as  the  one 
that  has  formed  a  covalent  bond  with  a  surface  (Ref  12  and  36).  For 
this  reason  bonding  and  chemisorption  are  hereafter  used  interchangeably. 
Wells  (Ref  30)  gives  the  S1-0  bond  energy  In  the  S102  molecule  as  108 
kcal . 
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Chemisorption  Data 

In  Chapter  II,  the  sticking  coefficient  was  expressed  as 

S  -  SQ  (1  -  ©)  (11) 

This  form  of  S  has  been  experimentally  verified  by  Antonlnl  (Ref  37)  for 
C02  up  to  180°C  and  for  02  on  "fresh"  silica,  and  by  Christmann,  et. 
al.,  (Ref  38)  for  H2  on  nickel.  However,  oxygen  atoms  have  not  been 
observed  to  chemisorb  to  natural  silicon  dioxide  surfaces  (Ref  27). 
Further  attempts  to  see  oxygen  atoms  chemisorbed  to  natural  silicon 
dioxide  surfaces  have  failed  (Ref  39). 

Antonlnl  (Ref  37)  observed  adsorption  sites  on  freshly  created 
glass  and  silica  surfaces.  These  sites  do  not  exist  on  glass  and  silica 
surfaces  In  their  natural  form.  As  a  result,  he  was  able  to  empirically 
determine  the  Initial  sticking  coefficient  for  oxygen  molecules  on  fresh 
silica  as  a  function  of  temperature  from  300°K  to  800°K.  His  results 
may  be  expressed  as 

SQ  =  0.0375  exp  (-0.0021  T)  (45) 

Hochstrasser  and  Antonlnl  (Ref  40)  found  that  there  are  5  x  10^ 

2 

bonds/cm  in  the  bulk  of  silicon  dioxide.  This  value  will  be  used  later 
as  the  estimate  of  the  number  of  adsorption  sites  on  the  silicon  dioxide 
surface. 


•!V, 


w 


V.  Evaluation  of  Initial  Model  and  Proposal  of  a  New 
Oxygen  Recombination  Model 


In  Chapter  III  the  approach  of  Jumper,  et.  al.,  (Ref  14)  was  followed 
to  derive  an  expression,  equation  (44),  for  the  recombination  coefficient 
when  both  oxygen  atoms  and  oxygen  molecules  are  assumed  to  chemisorb  to 
the  silicon  dioxide  surface.  This  model  may  be  Investigated  to  determine 
the  influence  of  changing  the  values  of  the  interrelated  parameters  In 
terms  of  attempting  to  match  the  character  of  the  data  described  In  the 
previous  chapter. 


Sensitivity  Analysis 

Following  Jumper's  approach,  the  activation  energy  for  the  oxygen 


recombination  on  silicon  dioxide  was  assumed  negligible  which  eliminated 


the  exp  (-  -py)  factors  from  equation  (44).  The  recombination  coefficient 
B 


was  then  calculated  as  a  function  of  temperature  from  300°K  to  2000°K 
using  realistic  values  for  the  Initial  sticking  coefficients  for  the 
oxygen  atoms  and  oxygen  molecules,  number  of  surface  sites,  sterlc  factor, 
and  the  thermal  desorption  energies  for  the  oxygen  atoms  and  oxygen 
molecules.  The  values  for  these  six  parameters  were  assumed  to  be  Inde¬ 
pendent  of  temperature.  Further,  the  initial  sticking  coefficient  for 
the  atoms  was  first  assumed  equal  to  the  Initial  sticking  coefficient 
for  the  molecules.  A  study  of  the  Influence  of  changing  the  values  of 
the  parameters  was  made  In  an  attempt  to  match  the  data  from  Greaves  and 
llnnett  (Ref  17)  and  also  to  keep  the  recombination  kinetics  first  order 


with  respect  to  the  atom  pressure  at  low  temperatures  and  transition  to 
second  order  only  at  some  elevated  temperature. 

-4 

This  study  led  to  a  recombination  coefficient  of  4.6  x  10  at 
300°K  If  the  values  of  the  parameters  given  in  Table  I  were  used. 
Further,  the  recombination  coefficient  experienced  a  maximum  value  of 
9.7  x  10'2  at  940°K  and  then  decreased  to  1.4  x  IQ"2  at  2000°K,  leading 
to  a  recombination  process  which  was  first  order  with  respect  to  the 
pressure  of  the  oxygen  atoms  at  300°K,  transitioning  to  second  order  at 
a  temperature  of  approximately  1300°K. 

Using  the  values  for  parameters  as  given  in  Table  I,  one  is  also 


able  to  Infer  the  physical  processes  at  work  In  this  model  as  a  function 
of  temperature.  At  300° K,  the  fourth  term  in  the  denominator  of 

•  •  C 

equation  (44),  NQ  SQ  Nq  P  exp  (-  -py),  Is  much  larger  than  the  other 

Oa  2  B 

terms  in  the  denominator.  Since  the  Initial  sticking  coefficients  for 
atomic  and  molecular  oxygen  are  assumed  equal,  the  recombination 
coefficient  may  be  expressed  as 


The  recombination  coefficient  thus  Increases  as  the  thermal  desorption 
of  the  molecules  Increases.  The  model  thus  predicts  that  the  fractional 
surface  concentration  of  the  oxygen  molecules  Is  essentially  unity  at 
300° K  and  decreases  as  the  fractional  surface  concentration  of  the  atoms 
increases.  Note  that  this  model  with  the  values  of  the  parameters  as 
given  in  Table  I  has  the  recombination  process  first  order  with  respect 
to  the  atomic  species  but  inversely  proportional  to  the  pressure  of  the 
oxygen  molecules. 


Above  800°K  the  first  term  dominates  the  denominator  of  equation 
(44)  so  that 


Y  -  2P 


(47) 


Thus,  at  this  temperature,  the  recombination  kinetics  remain  first  order 
with  respect  to  the  atomic  pressure  but  become  independent  of  the 
molecular  oxygen  pressure. 


TABLE  I.  Computational  Parameters  for  Oxygen  Atom 
Recombination  on  Silicon  Dioxide  Assuming  Atomic 
and  Molecular  Oxygen  Chemisorption. 


PARAMETER 

VALUE 

P,  STERIC  FACTOR 

0.1 

E,  ACTIVATION  ENERGY 

0 

Sn  ,  INITIAL  STICKING  COEFFICIENT 
°0  FOR  ATOMIC  OXYGEN 

0.3 

Sn  ,  INITIAL  STICKING  COEFFICIENT 
°02  FOR  MOLECULAR  OXYGEN 

0.3 

C  ,  SURFACE  SITES 

a 

3.0  x  106 

Dn,  THERMAL  DESORPTION  ENERGY  FOR 
u  ATOMIC  OXYGEN 

0.8  x  10"19  J/atom  = 
11.5  kcal/mole 

Dn  ,  THERMAL  DESORPTION  ENERGY 
u2  FOR  MOLECULAR  OXYGEN 

0.32  x  10“19  J/atom  = 
4.6  kcal/mole 

As  the  temperature  continues  to  increase,  the  recombination 
coefficient  experiences  a  maximum  at  940°K,  and  at  1300°K  the  third  term 
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In  the  denominator  begins  to  dominate.  As  a  result. 


and  the  recombination  process  becomes  second  order  with  respect  to  the 
atomic  species. 

Thus,  the  model,  adapted  wholesale  from  Reference  14,  Is  capable  of 

matching  the  character  of  the  available  catalytic-oxygen-recombination 

data  If  the  values  of  the  parameters  In  Table  I  are  used.  These  values 

are  almost  reasonable  except  for  the  extremely  low  value  for  the  number 

14 

of  adsorption  sites.  There  are  approximately  10  potential  surface 
sites  per  square  centimeter  of  surface  (c.f. ,  Chapter  IV).  To  have  only 
10®  sites  active  seems  too  low I  Further,  the  thermal  desorption  values 
for  atomic  and  molecular  oxygen  appear  somewhat  low  since  usual  values 
of  the  chemisorption  bond  energy  tend  to  be  between  15-200  kcal/mole 
(Ref  20). 

Finally,  the  assumption  that  oxygen  molecules  chemisorb  to  the 
silicon  dioxide  surface  Is  both  unlikely  from  the  data  (c.f..  Chapter 
IV),  and  not  in  keeping  with  the  usual  assumptions  concerning  molecular 
adsorption.  Laldler  (Ref  41),  for  example,  gives  a  convincing  argument 
which  concludes  that  molecular  chemisorption  In  general  is  negligible 
except  at  very  low  temperatures.  Ehrlich  (Ref  42)  also  explains  that 
molecular  chemisorption  is  not  observed  for  either  hydrogen  or  nitrogen 
on  many  metals. 


New  Model 


As  a  result  of  the  arguments  against  molecular  chemisorption,  the 
lack  of  experimental  evidence  that  molecular  oxygen  chemisorbs  to  natural 
silicon  dioxide  surfaces,  the  extremely  low  number  of  surface  sites 
required  for  this  model,  and  the  functional  dependence  of  the  recombina¬ 
tion  coefficient  upon  the  molecular  impingement  rate,  the  original  model 
was  set  aside.  A  new  model  was  developed  which  no  longer  required 
chemisorption  of  oxygen  molecules.  The  role  played  by  chemisorbed 
oxygen  molecules  was  replaced  by  hypothesizing  that  the  gas  phase  oxygen 
atoms  combine  directly  with  the  oxygen  atoms  that  constitute  the  silicon 
dioxide  surface.  Although  a  rather  radical  departure  from  previous 
approaches.  It  Is  really  only  a  slight  modification  of  thinking  in  terms 
of  the  Langmulr-Rldeal  recombination  mechanism  where  the  surface  atoms 
now  become  the  active  surface-atom  participants.  In  this  new  scheme 
then,  gas  phase  oxygen  atoms  replace  the  oxygen  atoms  as  they  are  removed 
from  the  surface-S1-0  lattice  by  recombination  reactions.  In  fact,  this 
idea,  although  not  used  extensively,  Is  not  without  precedent;  Llnnett 
and  Marsden  (Ref  27)  proposed  the  same  idea  in  1956  and,  as  we  will  see, 
the  experimental  evidence  seems  to  support  the  hypothesis. 

If  one  now  adapts  this  new  model,  a  pattern  begins  to  emerge  from 
the  data  which  at  least  qualitatively  supports  the  notion.  The  results 
of  Rosner's  work  (Ref  9),  for  example,  can  be  interpreted  as  support  for 
the  Idea  that  the  surface  sites  are  the  silicon  atoms  rather  than  the 
oxygen  atoms.  For  nitrogen  atom  recombination  on  silicon  dioxide,  Rosner 
found  nitrogen  atoms  Incorporated  into  the  S102  surface  as  silicon  nitrides 
or  silicon  oxynitrides.  This  experimental  data  support  the  hypothesis 


35 


that  the  silicon  atoms  on  a  silicon  dioxide  surface  are  the  active  surface 
sites.  Oxygen  atoms  bond  to  these  silicon  atoms  forming  the  sillcon- 
dloxlde-matrlx  surface  In  a  manner  not  unlike  chemisorption  of  oxygen  on 
a  silicon  surface.  Gas  phase  oxygen  atoms  could  then  be  colliding  and 
combining  with  the  oxygen  atoms  which  form  the  surface  of  the  silicon 
dioxide  Itself,  with  atoms  from  the  gas  phase  then  replacing  the  atoms 
lost  from  the  surface. 


In  the  new  model  the  oxygen  atom  that  Is  on  the  surface  and  possibly 
doubly  bonded  to  the  silicon  atom  (c.f. ,  Chapter  IV)  Is  the  oxygen  atom 
that  combines  via  the  Langmulr-Rideal  mechanism  with  a  gas  phase  oxygen 
atom.  The  silicon  atom  Is  the  surface  site  to  which  another  gas  phase 
oxygen  atom  readily  attaches  as  follows: 
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Rate  Equation 

The  time  rate  of  change  of  the  number  of  surface-matrix-bonded 
oxygen  atoms  per  unit  area  when  atoms  are  chemisorbing,  thermally  desorbing, 
and  recombining  via  the  Langmulr-Rideal  mechanism  and  subsequently  desorbing 


as  molecules,  can  now  be  expresses  as: 


=  N  SQ(1  -  0)  -  60  -  P  N  0  exp(-  /y)  (49) 

B 

In  a  steady  state,  this  equation  can  be  solved  for  the  fraction  of  the 
surface  that  Is  covered  by  oxygen  atoms  In  the  form: 

0  =  y  (50) 

S0N  +  5  +  PN  exp(-  -^j) 

B 


Recombination  Coefficient 

The  recombination  coefficient  is  then  given  by: 


Y  =  2P0  exp(-  yj) 
B 


If  we  substitute  equation  (50)  for  0, 


(51) 


2  P  SQ  N  exp(-  -jAr) 
y  ■ _ _ _ kb‘ 

SQ  N  +  6  +  P  N  exp(--j^T) 

B 


(52) 


A  quick  check  of  equation  (50)  for  the  fractional  surface  concentra¬ 
tion  shows  that  under  steady  state  conditions  the  fractional  surface 
concentration  would  be  equal  to  i  if  both  the  initial  sticking  coefficient 
and  the  sterlc  factor  are  equal  to  unity  and  if  thermal  desorption  of 
the  atoms  Is  negligible.  This  means  that  under  the  conditions  for  maximum 
recombination  only  half  the  sites  are  filled  at  any  Instant  In  time.  A 
steady  state  value  of  i  for  0  gives  a  maximum  value  of  unity  for  gamma, 
i.e.,  a  fully  catalytic  surface. 

A  qualitative  look  at  the  terms  In  the  denominator  of  equation  (52) 


for  the  recombination  coefficient  is  very  instructive  (Ref  18).  If  the 


term  SQN  is  large  in  comparison  to  the  other  two  terms,  then  equation 


(50)  reveals  that  0*1  and  that: 


y  *  2  P  exp(-  py) 
B 


(53) 


Since  the  recombination  coefficient  is  independent  of  the  impingement 
rate,  the  recombination  process  would  be  first  order  with  respect  to  the 


atom  pressure.  Also  note  that  the  recombination  process  is  controlled 


by  the  recombination  desorption  rate.  If  P  N  exp(-  py)  Is  large,  then  0 

B 

will  be  small  and: 


Y  =  2  S, 


(54) 


Then  the  recombination  process  would  again  be  first  order,  but  now  the 
recombination  process  is  controlled  by  the  chemisorption  rate.  If  5  is 
large,  then  0  will  be  small  and: 


2  P  50  N  exp{-  t^) 


(56) 


The  recombination  process  is  now  second  order. 

In  order  for  the  fractional  surface  coverage,  9,  to  be  large  at  low 


temperatures,  this  model  requires  SQ  N  to  be  large  at  these  low  tempera¬ 


tures;  no  other  option  is  possible.  However,  at  elevated  temperatures 
where  0  may  be  smaller,  the  order  of  the  recombination  process  must  be 
determined  empirically  in  order  to  know  which  term  will  dominate  the 
denominator  of  equations  (50)  and  (52).  If  the  order  is  first  at  elevated 
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temperatures,  then  P  N  exp(-  -py)  must  dominate.  For  second  order 

B 


recombination  at  the  higher  temperatures,  this  model  requires  high 
thermal  desorption  rate  for  the  atoms. 

Equation  (52)  for  the  recombination  coefficient  as  a  function  of 
temperature  requires  data  for  the  steric  factor,  sticking  coefficient, 
number  density,  apparent  activation  energy,  number  of  adsorption  sites, 
and  the  chemisorption  energy.  The  significance  of  the  model  Is  that  the 
recombination  coefficient  Is  not  dependent  upon  all  of  these  variables 
at  all  temperatures.  For  a  specific  species  of  gas  and  surface  at  low 
temperatures  where  0=1  and  the  recombination  process  Is  first  order, 
an  Arrhenius  plot  of  the  log  of  Y  versus  the  reciprocal  of  absolute 
temperature  will  provide  sufficient  data  so  that  the  steric  factor  and 
the  apparent  activation  energy  can  be  determined.  If  the  recombination 
process  remains  first  order  at  elevated  temperatures  where  0  is  small, 
the  initial  sticking  coefficient  is  the  only  dependent  variable  in  the 
expression  for  the  recombination  coefficient.  On  the  other  hand,  if  the 
recombination  process  goes  second  order  at  elevated  temperatures, 
reasonable  estimates  can  be  obtained  for  the  sticking  coefficient  and 
chemisorption  energy  when  recombination  coefficient  data  are  available. 


Independent  experiments  like  chemisorption  and  thermal  atom  desorption 
studies  would  confirm  the  values  obtained  for  the  dependent  variables 
from  the  recombination  coefficient  results.  Since  all  of  the  necessary 
data  are  not  available  to  date  to  completely  determine  all  the  physio- 
chemical  parameters,  it  is  only  possible  to  show  how  the  Langmulr-Rideal 
mechanism  satisfactorily  models  the  recombination  of  oxygen  or  silicon 
dioxide,  and  to  give  the  most  reasonable  estimate  of  the  parameters. 


Details  of  Model  for  Oxygen  Recombination  on  Silicon  Dioxide 

Since  in  this  model  the  surface  concentration  is  that  of  the  Si-0 
surface,  0  is  going  to  be  large  at  lower  temperatures  (0  -  1).  At  the 
low  temperatures,  then,  the  recombination  coefficient  is  given  by: 

Y  -  2  P  exp(-  -py)  (53) 

B 

If  both  the  steric  factor  and  the  activation  energy  are  independent  of 

temperature,  a  plot  of  the  log  of  the  recombination  coefficient  versus 

the  reciprocal  of  the  absolute  temperature  would  give  a  straight  line. 

However,  as  given  in  the  previous  chapter,  Llnnett's  data  increases  In 

value  nonli nearly  as  the  temperature  increases.  Greaves  and  Linnett 

(Ref  17)  proposed  that  the  activation  energy  increases  from  1  kcal/mole 

at  room  temperature  to  13  kcal/mole  at  300eC.  If  13  kcal/mole  is  inserted 

into  equation  (53)  along  with  a  value  of  1.6  x  10  for  the  recombination 

coefficient  at  300°C,  the  steric  factor  must  equal  73!  This  Is  not 

acceptable  since  by  definition  P  must  be  £  1.0.  If,  however,  the  steric 

factor  is  assumed  to  Increase  with  temperature,  and  the  activation  energy 

is  taken  to  be  constant  at  Its  room  temperature  value  of  1  kcal/mole 
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(6.949  x  10  J/atom),  a  very  satisfactory  match  of  Llnnett's  data  is 
obtained  for  an  exponentially  Increasing  steric  factor  with  temperature 
of  the  following  form: 

P  =  0.0000224  exp  (0.00908  T)  (56) 

This  value  increases  from  0.0C03414  at  300°K  to  0.1  at  925°K.  At  higher 
temperatures  the  steric  factor  is  held  equal  to  0.1. 


Assessing  the  Reasonableness  of  the  Functional  Dependence  of  the  Steric 
Factor  on  Temperature 


The  steric  factor  is  often  explained  as  the  fraction  of  the  suffi¬ 
ciently  energetic  collisions  that  have  the  proper  orientation  for  reaction 
to  occur  (Ref  43-44).  Adamson  states  that  the  orientation  explanation 
is  probably  valid  for  steric  factors  between  1.0  and  0.1  but  would  not 
likely  be  valid  for  very  small  values. 

The  term  exp  (-  -py)  represents  the  fraction  of  collisions  that 
B 

have  sufficient  energy  E  to  react.  This  assumes  that  E  comes  only  from 
the  kinetic  energy  of  relative  motion  of  the  colliding  particles. 

However,  If  the  internal  energy  of  the  particles  could  provide  the  activa¬ 
tion  energy  E  (Ref  44  and  22),  or  if  the  relative  translational  motion 
perpendicular  to  the  line  of  centers  could  provide  E,  then  we  must  modify 


exp(-  i~y)  as  follows: 
B 


Probability  that  f 
degrees  of  freedom 
have  sufficient 
energy  E 


ivi 


4 


i) 
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exp(-  Tc^T 


)  (571 


The  steric  factor  P  could  then  be  defined  to  include  this  pre-exponential 
factor: 

E  1  (?  ‘  11 

s! _ 

4  -  1)! 


which  Is  a  function  of  both  temperature  and  the  number  of  degrees  of 
freedom  available  to  supply  E  (Ref  45).  Vlncentl  and  Kruger  explain 
that  only  a  rough  estimate  can  be  made  of  f,  due  to  the  doubt  as  to 


whether  the  vibration  of  the  catalytic  molecule  can  be  effective  and 
whether  the  relative  translational  motion  normal  to  the  line  of  centers 
should  be  considered.  Although  doubts  exist,  the  sterlc  factor  could  be 
temperature  dependent. 

Another  explanation  of  the  sterlc  factor  Is  given  by  Clarke  and 
McChesney  in  terms  of  partition  functions  of  the  reactants: 

TWO  ATOMS  P  =  1 

Atom  +  Diatomic  molecule  Nonlinear  complex  P  =  (Z^l (zrot)_1 

Linear  complex  P  =  (zv1b)2(zrot)"2 

3  3 

Two  diatomic  molecules  Nonlinear  complex  P  =  (zv-fb)  (zrot) 

Linear  complex  P  =  (zvib)4(zrot}~4 

They  explain  that  at  ordinary  temperatures  Zv^b  *  1  and  ZrQt  <  so 
that  for  reactions  between  two  atoms  or  an  atom  and  a  diatomic  molecule, 
the  steric  factor  is  between  0.1  and  1.  As  the  reacting  molecules  become 
more  and  more  complicated,  the  sterlc  factor  decreases  In  value. 

Even  though  these  explanations  give  a  functional  dependence  of  the 
sterlc  factor  on  temperature,  the  actual  relationship  depends  on  the 
assumptions  made  and  can  only  be  verified  by  further  detailed  experiments. 
Therefore,  with  the  data  that  is  available,  the  explanation  given  by 
Greaves  and  Llnnett  (Ref  17)  is  still  probably  the  best.  They  state 
that  the  nonlinear  Increase  In  the  recombination  coefficient  with  temp¬ 
erature  (therefore,  nonlinear  increase  In  the  sterlc  factor  for  a  constant 
activation  energy)  could  be  explained  In  terms  of  changing  the  number  of 
active  surface  sites.  Dickens  and  Sutcliffe  (Ref  13)  found  that  the 


recombination  coefficient  increased  nonlinearly  with  temperature  for 


MngOg*  FegO^,  N10,  CuO,  CdO,  and  C^O^.  They  explained  this  behavior  as 
being  due  to  a  change  in  state  of  the  surface.  The  relationship  among 


the  crystalline  forms  of  silicon  dioxide  is  as  follows  (Ref  46): 


«  quartz 


573°C  870°C  1470°C  1710°C 

- — -  8  quartz  - — -  8  tridymlte  - — *•  6  Cristobal ite  — *  liquid 

4 1  160  to  "*  A  |  200  to  "r 

I  i  120°C  I  I r  275°C 

“  tridymlte  *  cristobal ite 


Since  silicon  dioxide  changes  its  surface  structure  with  temperature,  it 


appears  that  the  best  explanation  for  the  steric  factor  increasing  with 


increasing  temperature  is  that  the  number  of  active  surface  sites  Increases 


with  temperature  as  the  surface  structure  changes  causing  the  steric 


factor  to  also  Increase. 


In  the  present  model  the  surface  sites  are  taken  to  be  constant 


with  temperature.  Thus,  if  the  explanation  of  Greaves  and  Linnett  is 


correct,  one  must  interpret  the  steric  factor  to  include  the  surface 


site  change  as  well  as  its  normal  meaning.  It  must  be  recognized,  then. 


that  the  definition  used  here  is  at  slight  variance  with  the  classical 


interpretation.  What  has  not  changed,  however,  is  the  real  fact  that 


the  steric  factor  has  always  been  and  continues  to  be  here  a  catch-all 


factor  which  includes  the  uncertainties  In  the  model,  the  point  being 
made  that  it  is  not  unreasonable  that  it  has  a  temperature  dependence. 


Further  Details  of  Model 


Since  the  recombination  process  remains  first  order  to  at  least 


1100°K  which  is  above  the  temperature  where  gamma  is  a  maximum,  the 


recombination  coefficient  is  approximated  by  the  following: 


Y  -  2  SQ  (54) 

Figures  3  and  4  show  how  the  Initial  sticking  coefficient  Is  only  of 
significance  at  elevated  temperatures.  The  influence  of  the  quantitative 
changes  In  the  Initial  sticking  coefficient  on  the  recombination  coeffi¬ 
cient  is  also  shown.  The  same  functional  form  as  Antoninl's  results 
(c.f..  Chapter  IV)  for  02  chemisorption  on  "fresh"  silicon  dioxide  is 
adapted  here  and  is  as  follows: 

SQ  =  0.05  exp (-0.002  T)  (58) 

for  the  Initial  sticking  coefficient  of  oxygen  atoms  onto  silicon 
dioxide.  The  Initial  sticking  coefficient  decreases  from  0.02744  at 
300° K  to  0.0009158  at  2000°K. 

At  a  temperature  somewhere  above  1100°K,  the  thermal  desorption  of 
atomic  oxygen  becomes  significant.  Equation  (13)  gives  the  thermal 
desorption  rate  per  unit  area  of  covered  surface  as: 

C  kRT  D 

{ •  -V-  «p<-  tjr>  (13> 

C  Is  the  number  of  adsorption  surface  sites  per  unit  area,  and  Hochstrasser 

a 

and  Antonlnl  (Ref  40)  found  that  there  are  5  x  10*4  bonds/cm2  In  the 
bulk  of  silicon  dioxide.  Therefore,  I  will  use  this  value  as  the  number 
of  adsorption  sites  on  the  surface  of  silica. 

Wells  (Ref  30)  gives  the  Si-0  bond  energy  In  the  SIO^  molecule  as 
108  kcal.  Therefore,  each  half-bond  strength  would  be  54  kcal.  As  a 
reasonable  estimate,  150%  of  the  half-bond  strength  is  used  as  the  energy 


required  to  extract  from  the  surface  the  oxygen  atom  that  has  the  double 
bond  to  the  silicon  atom.  This  81  kcal  (5.629  x  10"^  J/atom)  Is  the 
value  used  for  the  chemisorption  well  depth.  Therefore,  all  variables 
are  quantified  in  equation  (52)  for  the  recombination  coefficient  as  a 
function  of  temperature.  Figure  5  demonstrates  the  Influence  of  the 
thermal  desorption  energy  on  the  recombination  coefficient.  The  desorp¬ 
tion  energy  was  decreased  in  5  kcal/mole  Increments  from  85  to  45  kcal/mol 

For  negligible  second  order  recombination  at  1100°K  the  thermal  desorption 
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energy  must  be  at  least  60  kcal/mole  (4.169  x  10  J/atom).  An  experi¬ 
mental  investigation  as  to  the  temperature  where  the  recombination  process 
turns  to  second  order  would  be  very  useful. 

Figure  6  gives  the  results  of  the  model  for  the  recombination  of 
oxygen  atoms  on  silicon  dioxide  using  the  values  for  the  parameters 
summarized  In  Table  II.  A  close  look  at  the  calculations  performed  In 
20°K  Increments  from  300°K  to  2000°K  is  very  informative.  The  physical 
and  chemical  microscopic  details  are  readily  available  from  this  model. 

As  was  discussed  earlier,  the  S0N  term  In  the  denominator  of  equation 
(52)  for  the  recombination  coefficient  must  be  the  dominant  term  at  room 
temperature  in  order  for  the  recombination  process  to  be  first  order  and 
for  the  fractional  surface  concentration  to  be  approximately  unity. 

From  the  details  of  the  model  this  term  dominates  from  300° K  to  740° K. 


With  the  SqN  dominating  In  this  temperature  range,  the  recombination 
process  is  controlled  by  the  recombination  desorption  rate  through  the 
steric  factor  and  the  activation  energy.  From  740°K  to  1600°K,  the 


recombination  process  remains  first  order,  and  the  recombination  coeffi- 

clent  has  a  maximum  value  of  1.4  x  10  at  900°K.  The  P  N  exp  (-  -j^y) 

B 
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term  In  equation  (52)  dominates  In  this  temperature  range  with  the  recom¬ 
bination  process  being  controlled  by  the  chemisorption  rate  through  the 
Initial  sticking  coefficient.  Above  1600°K  the  thermal  desorption  rate 
term  Is  the  largest  In  the  denominator  of  equation  (52).  At  these  temp¬ 
eratures  the  recombination  process  Is  second  order,  and  the  fractional 
surface  concentration  Is  very  small.  The  oxygen  atoms  that  are  surface- 
matrix-bonded  rapidly  leave  the  surface  again  due  to  the  high  thermal 
desorption  rate.  There  simply  aren't  many  surface  sites  filled  with  a 
surface-matrix-bonded  oxygen  atom,  which  Is  a  prerequisite  for  heterogeneous 
recombination.  As  a  result,  the  recombination  coefficient  is  very  small 
and  continues  to  decrease  In  value  as  the  temperature  Increases. 


TABLE  II.  Computational  Parameters  for  Oxygen  Atom  Recombination 
on  Silicon  Dioxide  Assuming  Atomic  Oxygen  Chemisorption  Only. 


PARAMETER 
P,  STERIC  FACTOR 

E,  ACTIVATION  ENERGY 

SQ,  INITIAL  STICKING  COEFFICIENT 
C  .  SURFACE  SITES 

a 

D,  THERMAL  DESORPTION  ENERGY 


VALUE 

0.0000224  exp  (0.00908  T) 
(maximum  value  *  0.1) 

1  kcal/mole  = 

6.949  x  10" J/atom 

0.05  exp  (-0.002  T) 

5  x  1014  sltes/cm^ 

81  kcal/mole  =,Q 

5.629  x  10"*y  J/atom 


The  fractional  surface  concentration  Is 
function  of  the  temperature.  The  adsorption 
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shown  In  Figure  7  as  a 
rate,  thermal  desorption 


*  JK  >. 

it 


rate,  and  the  rate  of  recombination  desorption  are  shown  In  Figure  8. 

The  total  rate  of  removal  of  chemisorbed  (surface-matrix-bonded)  atoms 
Is  the  sum  of  the  thermal  desorption  rate  and  the  recombination  desorption 
rate.  At  steady  state  conditions  this  total  removal  rate  equals  the 
adsorption  rate. 


VI.  Model  Demonstration 


In  order  to  demonstrate  the  use  of  the  present  recombination  model 
and  to  explore  the  rational  limits  which  may  be  placed  on  the  role  of 
catalysis  for  an  engineering  problem,  a  set  of  boundary-layer  calculations 
were  made  using  numerical  methods.  As  indicated  in  the  introduction, 
that  portion  of  the  Space  Shuttle  reentry  where  oxygen  Is  dissociated 
but  nitrogen  Is  not  was  chosen  for  an  example  problem.  Such  a  portion 
of  the  Shuttle  reentry  is  below  200,000  feet.  According  to  Dunn  and 
Kang  (Ref  47),  this  allows  the  analysis  of  a  one-foot  nose  radius  body 
(of  which  the  Shuttle  is  one)  to  be  treated  by  conventional  boundary- 
layer  analysis.  Above  200,000  feet  an  alternate  method  Is  required  due 
to  the  thickening  of  the  viscous  region  within  the  shock  layer. 

To  suggest  that  the  analysis  described  below  exactly  models  the 
Space  Shuttle  problem  is  misleading,  but  a  geometry  somewhat  akin  to 
the  Shuttle  geometry  was  analyzed.  This  geometry  was  that  of  a  1.0 
foot  radius  spherical  Reinforced  Carbon-Carbon  (RCC)  nose  blended  to  a 
40°  half-angle  cone.  This  cone  Is  made  of  High-Temperature  Reusable 
Surface  Insulation  (HRSI)  tiles  which  are  coated  with  silicon  dioxide. 
Figure  9  is  a  schematic  of  the  modeled  geometry  indicating  the  compu¬ 
tational  domain.  This  domain  is  from  the  X  =  61.  cm.  location  to 
approximately  the  x  =  90.  cm.  location.  This  chapter  contains  an 
explanation  of  the  numerical  boundary-layer  algorithm  that  was  used  and 
the  results  obtained. 


Rei nforced 
Carbon-Carbon  (RCC) 


Boundary-Layer  Code 

After  the  61.  cm.  point,  the  boundary-layer  was  modeled  using  a 
two-dimensional,  rectilinear  computational  space,  but  undergoing  a  pressure 
distribution  for  the  axisymmetrlc  cone  (c.f. ,  above  and  below).  The 
boundary-layer  code  is  a  conventional  two-dimensional  boundary-layer 
method  which  has  been  described  elsewhere  (Ref  8),  where  it  was  validated 
for  a  high  temperature,  compressible,  reacting  flow  (including  catalytic 
walls)  of  fluorine  atoms,  molecules  and  diluent.  The  code  was  further 
validated  in  Reference  48  for  low  temperature,  compressible  air.  The 
equations  for  continuity,  momentum,  species,  and  energy  in  boundary-layer 
form  from  Reference  8  are  as  follows: 


Continuity:  -g|(pu)  +  -gp(pv)  = 

0 

(59) 

Momentum:  pu|£  +  pv 

3yv  3y' 

(60) 

3w.  3w. 

Species:  pih^  +  Pv  gy  =  mj‘ 

a 

' '  +  W(p  DS^yL) 

(61) 

Energy:  0Cp  («£♦<$-«£♦  „<# •)*  -  >  - 

3w.  3h.  . 

2  p  DJlyV  ’ 2  hJ  V  "  (62) 


where  u  is  the  x  component  of  velocity,  v  is  the  y  component  of  velocity,  p 
is  the  mixture  density,  p  is  the  pressure,  y  is  the  mixture  viscosity, 

Wj  is  the  mass  fraction  of  species  j,  m j ’ * '  is  the  rate  of  production  of 
species  j  per  unit  volume,  D.  is  the  diffusivity  of  species  j  into  the 


mixture,  Cp  Is  the  specific  heat  at  constant  pressure  for  the  mixture,  T 


is  the  temperature,  h.  is  the  absolute  enthalpy  for  species  j,  and  k  is 

J 


the  thermal  conductivity  for  the  mixture.  The  von  Mises  transformation 


was  then  used  where  \[>  is  the  compressible  stream  function  given  by: 


M.  pu 
3y  M 


dip 

Tx  -pv 


The  transformed  equations  are  then  given  by  the  following: 


Momentum:  -g£  -  -  -ft  +  -gffrPU  |^) 


Dj  p2  “3F  +  ^(kpu^) 


An  explicit  finite  difference  method  as  given  in  Reference  8,  was  used 


to  solve  equations  (65)  -  (67). 


Thermodynamic  and  Transport  Properties 


The  sensible  enthalpy  and  specific  heats  for  each  of  the  five  species 
were  extracted  from  the  JANAF  data  (Ref  49).  The  chemical  enthalpy  of 
formation  from  the  JANAF  data  was  added  to  the  sensible  enthalpy  to 


obtain  the  absolute  enthalpy.  The  specific  heat  for  the  mixture  is 


^ i) >)Y 'i i  l \ 


calculated  as  follows: 


N 

Cn  =  Z  x<  Cn 
p  j-i  3  pj 


(68) 


where  x.  is  the  mole  fraction  of  species  j  and  C  is  the  specific  heat 

J  PJ 

of  species  j.  The  mole  fractions  were  calculated  from: 

w. 


'j 


(69) 


where  M.  and  M  are  the  molecular  weights  of  the  individual  species  j  and 

J 

mixture,  respectively,  and  w.  Is  the  mass  fraction  of  species  j.  The 

J 

molecular  weight  M  is  calculated  as  follows: 


The  mixture  density  was  calculated  by: 

(7D 

where  p  Is  the  absolute  pressure,  T  is  the  absolute  temperature  and  R  is 
the  universal  gas  constant.  The  viscosity,  thermal  conductivity,  bimo- 
lecular  dlffusivlty,  and  the  diffusivity  of  species  1  into  the  mixture 
were  calculated  as  explained  by  Jumper,  et.  al.,  (Ref  8)  and  are  summarized 
here.  The  viscosity  was  calculated  as  (Ref  50): 

..  2.6693  x  urf/fiT 

y-j  - _ y  i _  (72) 

2 

ai  % 

where  ftp  Is  the  reduced  collision  integral  curve  fit  as  a  function  of 
kBT/er  The  va1ues  used  for  the  Lennard-Jones  (6-12)  characteristic 

energy  of  interaction  divided  by  the  Boltzmann  constant,  and  c^,  the 
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Lennard-Jones  (6-12)  characteristic  diameter,  came  from  Reference  51  and 
are  given  in  Table  III.  The  viscosity  of  the  mixture  was  calculated 


i 


CM 


using: 


^mix  ~  2 
i=l 


j=l  XJ  *1j 


,  M.  -i  I*,  1  M.  i  2 

*—  1  + 1 r  1  +tt  r 

1J  /e  "j  wj  Mi 


The  thermal  conductivity  was  calculated  using  the  following  relationships 
from  Reference  50. 


ki  =  (cp<  +i{r}  y1 


k  .  =  l  xik1 

mix  i=l  “N 

E  x, <!>,,• 


The  following  relationship  from  Reference  50  was  used  to  calculate  the 
blmolecular  diffuslvlty: 


0.0018583  lA^TT  +  lb) 


0°i }  nD 


where  was  curve  fit  as  a  function  of  kgT/e^  where 

*  (ei£/ 


and 


ou  *  (oi  +oj)/2 


(79) 


The  diffusivity  of  species  1  into  the  mixture  was  calculated  with 
(Ref  50): 


1  -  x 


D<  - 


T 

E 

j=l  D 


I 


(80) 


ij 


TABLE  III.  Lennard-Jones  (6-12)  Data 


Species 

_fl_ 

kB 

00 

a1 

(Angstroms) 

N 

71.4 

3.298 

0 

106.7 

3.050 

n2 

71.4 

3.798 

°2 

106.7 

3.467 

NO 

116.7 

3.492 

Homogeneous  Chemistry 

The  gas-phase  (homogeneous)  chemical 

reactions  considered 

research  are  shown  in  Table  IV.  These  are  the  same  chemical  reaction 
data  used  by  Kang  and  Dunn,  and  Miner  and  Lewis  for  a  five  species  gas 
(Ref  47,  52-53).  Note:  T"**^  was  dropped  from  Kang  and  Dunn's  equation 
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m 


6  In  Table  I  of  Reference  52.  The  VSLNQH  code  (c.f. ,  below),  used  to 
obtain  the  Initial  conditions,  uses  chemical  reaction  data  of  a  different 
form.  The  significance  of  this  difference  was  not  investigated. 


These  reactions  were  then  used  to  compute  the  rate  of  creation/deple 
tlon  in  the  following  manner  (Ref  22): 

d[N]  =  2kf  [N2]  [N]  -  2kb  [N]2  [N] 
at  6  6 

+2kf^  [N2]  [0]  -  2k^  [N]2  [0] 

+2kf  [N2]  [N2]  -  2kb  [N]2  [N2] 

8  8 

+2kf^  [N2]  [02]  -  2kbg  [N]2  [02] 

+2kf  [N?]  [NO]  -  2kh  [N]2  [NO] 

T10  2  D10 

+kf  [NO]  [N]  -  kh  [N]  [0]  [N] 

T11  D11 

+kf  [NO]  [0]  -  kh  [N]  [0]  [0] 

+kf  [NO]  [N?]  -  k.  [N]  [0]  [Np] 

T13  2  D13  2 

+kf  [NO]  [0?]  -  k.  [N]  [0]  [02] 
f14  2  bl4  2 

+kf  [NO]  [NO]  -  kh  [N]  [0]  [NO] 

T15  d15 

+kf  [NO]  [0]  -  kh  [Op]  [N] 

T16  D16  2 

+kf  [Np]  [0]  -  k.  [NO]  [N]  (81) 

t17  2  °17 
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d[N9] 


2\ 

C02]  [N]  - 

2kh 

bl 

[O]2  [N] 

+2k. 

t2 

[02]  [0]  - 

2kh 

b2 

to]2  [0] 

+2k- 

t3 

[02]  [N2] 

-  2kb 

[O]2  CM,] 

3  2 

+2k, 

t4 

[02]  [02] 

-2kb 

[O]2  [02] 

4  2 

+2k- 

t5 

[02]  [NO] 

‘2kbl 

[O]2  [NO] 

5 

+kf 

[NO]  [N]  - 

kb 

[N]  [0]  [N] 

T11 

bll 

+kf 

[NO]  [0]  - 

kb 

[N]  [0]  [0] 

'12 

b12 

+kf 

t13 

[NO]  [N2] 

-  K  [N]  [0]  [N?] 
°13  2 

+kf 

T14 

[NO]  [02] 

-  k.  [N]  [0]  [0?] 
D14  2 

+kf 

[NO]  [NO]  ■ 

“  kK 

CM3  [0]  [NO] 

t15 

15 

-k 

[NO]  [0]  + 

kb 

[02]  [N] 

16 

b16 

"kf17 

[n2]  [0]  + 

kh 

b17 

[NO]  [N] 

"kf6 

[n2]  [N]  + 

kb  [N]2  [N] 

°6 

"kf7 

[n2]  [0]  + 

kh 

b7 

[N]2  [0] 

“kf 

t8 

[N2]  [N2]  h 

h  kh 
b8 

[N]2  [Mj] 

~kf 

t9 

[N2]  [02]  h 

h  kh 

b9 

[N]2  [02] 

"kf 

T10 

[N2]  [NO] 

+  kb  [N]2  [NO] 
D10 

"kf 

t17 

[n2  [0]  + 

K 

b17 

[NO]  [N] 

L  2J  .  -kf  [Op]  [N]  +  kh  [C]"  [N] 
~St~  "  T1  L  D1 

-kf  [Op]  [0]  +  kb  [O]2  [0] 
f2  L  b2 

-kf  [o2]  [n2]  +  kb  [O]2  [n2] 

3  3 

-kf  C02]  [02]  +  kb  [0]2  [02] 

4  4 

-kf  [02]  [NO]  +  kb  [0]2  [NO] 

5  5 


+kf 

'16 

[NO] 

[0]  - 

kb  C°2] 
d16  c 

[N] 

-  -kf 

[NO] 

[N]  + 

kh  [N]  | 

:0]  [N] 

fll 

bll 

-kf 

[NO] 

[0]  + 

kb  [N]  | 

[0]  [0] 

f12 

b12 

”kf 
’  13 

[NO] 

[N2] 

+  kb  [N] 
b13 

[o]  [n2] 

-kf 

f14 

[NO] 

[02] 

+  kb  [N] 
b14 

[0]  [o2] 

-kf 

[NO] 

[NO] 

+  kb  [N] 

[0]  [NO] 

'15 

b15 

-kf 

[NO] 

[0]  + 

kb  ^2^ 

[N] 

T16 

b16  c 

+kf 

[n2] 

[0]  - 

kb  [NO] 

[N] 

f17 

b17 

(84) 


(85) 


where  the  subscript  on  the  forward  rate  constant  k^  and  backward  rate 
constant  kb  refers  to  the  reaction  In  Table  IV.  Is  the  rate  of 

3  •  i » i 

production  of  species  j  in  moles/cm  -sec,  m^  ,  that  Is  found  in  both 
the  species  and  energy  equations,  [j]  is  the  molar  concentration  of 

3 

species  j  in  moles/cm  and  is  calculated  as  follows: 


8i 


Boundary  Conditions  Due  to  Recombination 

At  the  wall,  the  rate  of  diffusion  of  atomic  oxygen  toward  the  wall 
must  just  equal  the  rate  at  which  atomic  oxygen  disappears  due  to  catalytic 
recombination  forming  molecular  oxygen.  Also,  the  rate  at  which  molecular 
oxygen  is  formed  determines  the  rate  at  which  the  molecular  oxygen  diffuses 
away  from  the  wall  (Ref  8).  This  results  in  the  following  boundary 
conditions: 


w0U)  = 


°0  w0  ^ 
kDt 


(87) 


Dn  +  Y  (  *B'  A 


2  urn, 


■)*  Ay 


wn  (!)  =  Wn  (2)  + 


wn  (1)  kRT  , 

-V 


(88) 


where  Wq(1)  and  Wq(2)  are  the  mass  fractions  for  atomic  oxygen  at  the 

wall  and  at  a  distance  Ay  from  the  wall,  respectively,  and  wn  (1)  and 

u2 

wn  (2)  are  the  mass  fractions  for  molecular  oxygen  at  the  wall  and  at 
u2 

the  first  cell  out  from  the  wall. 

As  discussed  by  Jumper,  et.  al.,  (Ref  8),  the  heat  transfer  enhance¬ 
ment  per  unit  area  due  to  the  recombination  of  oxygen  atoms  may  be 
expressed  as  the  following: 

qrecomb  P  Y  w0^  ^2  ^m^  (ho  "  h02^ 


where  the  enthalpies,  hg 


and  hn  , 
u2 


are  on  a  per  mass  basis. 


CC 


L*l«] 


The  surface  temperature  was  treated  in  two  ways:  the  wall  temperature 
could  be  specified  or  it  could  be  calculated  based  on  the  radiation 
equilibrium  assumption  (Ref  7).  Since  the  thermal  conductivity  of  the 
Space  Shuttle  Insulation  tiles  is  very  low,  l.e.,  only  a  small  amount  of 
energy  Is  conducted  through  the  surface,  the  surface  is  at  a  temperature 
high  enough  so  that  the  total  heating  is  assumed  equal  to  the  energy 
that  can  be  radiated  away  with  an  emissivity,  e  of  0.85: 

“total  ’  6  0  t4  <9°) 

where  q  Is  the  total  heating,  a  is  the  Stefan-Boltzmann  constant,  and  T 
is  the  surface  temperature.  When  the  second  option  was  used,  the  surface 
temperature  was  determined  by  the  total  heating  at  the  previous  step. 

Initial  Conditions 

The  code  written  by  Song  and  Lewis  (Ref  54),  VSLNQH:  An  Axisymmetric 
Viscous  Shock-Layer  Code  for  Noneciuil ibrium  Finite-Rate  Chemically  Reacting 
Viscous  Flows  over  Hyperboloid  Geometries,  was  first  used  to  obtain  the 
initial  conditions  for  the  boundary-layer  code.  The  VSLNQH  code  yielded 
mass  fractions,  temperature,  and  velocity  profiles  over  the  Reinforced 
Carbon-Carbon  (RCC)  nose  cap  to  approximate  a  Shuttle  flight  condition 
of  Mach  14.1  at  an  altitude  of  188,000  feet  assuming  the  RCC  surface  is 
at  a  uniform  temperature  of  1703°K  (3065CR)  and  is  noncatalytic.  (This 
VSLNQH  code  can  only  treat  noncatalytic  or  fully  catalytic  surfaces.) 

Table  V  and  Figures  10-12  show  the  Initial  temperature,  velocity,  and 
mass  fraction  profiles,  respectively,  at  the  streamwise  location  of  61.0 
cm.  This  corresponds  to  the  RCC  Interface  with  the  High-Temperature 
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Figure  10.  Initial  Temperature  Profile  at  x  =  61.0  cm 
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Figure  12.  Initial  Mass  Fraction  Profiles  at  x  =  61.0  cm. 


Reusable  Surface  Insulation  (HRS I )  tiles.  These  tiles  are  coated  with 
predominately  silicon  dioxide.  Figure  12  readily  shows  the  high  oxygen 
atom  content  (mass  fraction  of  19%)  of  the  flow  adjacent  to  the  surface 
due  to  this  flow  having  come  through  the  near-normal  portion  of  the  bow 
shock  and  also  not  having  been  reduced  by  the  noncatalytlc  assumption 
for  the  RCC  surface.  At  this  flight  condition  the  nitrogen  atom  content 
along  the  surface  is  negligible  (mass  fraction  of  N  =  0.0022%  at  x  = 

61.0  cm)  and,  as  such,  modeling  this  condition  allowed  for  independent 
consideration  of  the  oxygen  atom  recombination.  The  VSLNQH  code  is 
based  upon  the  assumption  that  the  Lewis  number  for  each  of  the  species 
is  equal  to  1.4: 


Ip  p  Cn 

j  - £ — *1—  assumed  ^  4 


(91) 


That  is,  once  the  density,  specific  heat,  and  thermal  conductivity  for 
the  air  mixture  is  calculated,  the  diffusivity  of  each  species  j  is 
calculated  using  Le.  =  1.4.  The  boundary-layer  code  used  in  this 

J 

research  does  not  make  this  assumption. 

The  VSLNQH  code  was  also  used  to  produce  an  initial  pressure  at  the 
RCC/HRST  interface  of  0.0424733  atm.  and  a  pressure  gradient  of: 

-t-  =  (3.74454  x  10"6)  x  -(4.31620  x  10‘4)  atm/cm  (92) 


valid  over  the  next  27  cm. 

Results  and  Discussion 

The  Song  and  Lewis  VSLNQH  code  was  run  with  the  RCC  surface  assumed 
to  be  noncatalytlc  and  at  1703°K.  This  code  was  further  run  to  x  *  89. 


cm  assuming  a  noncatalytic  surface  of  fixed  temperature.  The  surface 
temperature  was  set  at  1500°K  at  65.0  cm  and  1353°K  at  82.1  cm.  Linear 
temperature  distributions  were  assumed  between  these  temperatures  with 
the  temperature  remaining  1353°K  downstream  from  82.1  cm.  The  boundary- 
layer  code  was  then  run  with  the  initial  conditions  from  the  VSLNQH  code 
starting  at  the  RCC/HRSI  Interface  (x  =  61.0  cm)  with  this  same  surface 
temperature  distribution  and  with  y  =  0.00.  The  results  are  shown  in 
Figure  13.  The  results  for  the  explicit  boundary-layer  code  (y  =  0)  are 
significantly  higher  than  the  VSLNQH  code  results  (with  y  =  0)  Initially, 
but  at  89.0  cm  the  boundary-layer  code  results  are  7%  lower  than  the 
VSLNQH  code  results.  In  this  nose  region  heat  transfer  predictions 
often  differ  by  up  to  40%  (Ref  7).  For  comparison  the  results  of  the 
boundary-layer  code  with  Y  from  the  model  and  for  the  fully  catalytic 
surface  (y  =  1.00)  are  also  shown  in  Figure  13,  where  all  results  are 
for  the  specified  surface  temperature  distribution.  This  analysis  also 
assumes  that  all  of  the  energy  released  when  recombination  occurs  is 
transferred  to  the  surface. 

Heating  rates  are  shown  In  Figure  14  when  radiation  equilibrium  is 
assumed  on  the  HRS  I  tiles  for  three  runs  of  the  boundary- layer  code. 

The  corresponding  surface  temperatures  are  shown  in  Figure  15.  The 
STS-2  flight  test  datum  Is  shown  on  these  two  figures  only  for  comparison 
and  is  not  to  be  used  to  draw  a  conclusion  relative  to  the  catalycity  of 
the  surface. 

A  convergence  check  was  performed  by  approximately  halving  A\f>  and 
this  check  gave  an  error  of  0.9%  for  the  radiation  equilibrium  temperature 
at  80.  cm. 
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Figure  13.  Heating  Rate  for  Specified  Surface 
Temperature  Distribution. 


mva msm 


Figure  15.  Radiation  Equilibrium  Surface  Temperature. 
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A  significant  conclusion  that  can  be  drawn  from  Figure  14  Is  that 
the  heating  due  to  recombination  of  oxygen  atoms  can  be  large.  In 
comparison  with  the  heating  for  y  =  0,  the  Si  0^  surface  experiences  21% 
more  heating  and  a  fully  catalytic  surface  would  experience  33%  more  at 
89.  cm  when  radiation  equilibrium  is  assumed.  Even  though  Y  is  relatively 
small,  ranging  from  5.5  x  10  at  65.  cm  (1431°K)  to  7.2  x  10  at  89. 
cm  (1289°K),  the  total  heating  of  a  silicon  dioxide  surface  at  89.0  cm 
was  only  9.1%  below  what  would  have  been  experienced  by  a  fully  catalytic 
surface  at  these  simulated  flight  conditions.  Or,  put  another  way,  the 
heating  due  only  to  recombining  oxygen  predicted  by  the  model  is  approxi¬ 
mately  63%  of  that  which  would  be  predicted  for  y  =  1.0,  even  though  the 
Y's  differ  by  nearly  three  orders  of  magnitude.  This  can  be  explained 
by  noting  that  the  reaction  begins  to  become  diffusion  controlled  (i.e., 
a  limit  to  the  rate  at  which  oxygen  atoms  can  diffuse  to  the  surface) 
for  the  y  *  1.0  case. 

The  temperature  profiles  at  x  *  89.  cm  are  shown  in  Figure  16  for 
the  three  cases:  Y  =  0,  y  from  model,  and  y  *  1  when  radiation  equilib¬ 
rium  is  assumed  in  each  case.  No  significant  difference  exists. 

The  velocity  profiles  at  x  =  89.  cm  for  the  same  three  cases,  again 
with  the  radiation  equilibrium  assumption,  are  shown  in  Figure  17. 

There  is  no  significant  difference  in  these  profiles  either. 

The  initial  mass  fraction  profile  for  the  four  species  of  signifi¬ 
cance,  which  is  plotted  in  Figure  12,  has  ^W1  gradients  at  the  surface 

W 

equal  to  zero  for  all  the  species  since  this  profile  was  obtained  from 
the  VSLNQH  code  assuming  a  noncatalytlc  surface.  The  mass  fraction 
profiles  with  Y  from  the  model  at  the  downstream  location  is  shown  in 


Figure  18  where  the  oxygen  atom  and  molecule  mass  fractions  are  observed 
to  have  changed  significantly.  Note  also  the  increase  in  the  oxygen 
atom  concentration  out  in  the  free  stream  due  to  the  dissociation  of 
oxygen  molecules.  Figure  19  shows  the  oxygen  atom  and  molecule  mass 
fractions  for  the  three  cases.  The  only  significant  difference  in 
comparing  the  fully  catalytic  and  noncatalytic  mass  fraction  profiles  28 
cm  downstream  is  within  the  first  centimeter  from  the  surface.  The 
atomic  oxygen  mass  fraction  Is  not  zero  (0.07%  at  x  =  89.  cm)  at  the 
wall  even  when  the  surface  Is  modeled  as  fully  catalytic,  supporting  the 
comments  made  by  Jumper,  et.  al.,  (Ref  8). 

The  secondary  objective  of  this  research  has  been  to  demonstrate 
the  use  of  the  oxygen  recombination  model  with  proper  treatment  of 
surface  boundary  conditions.  This  demonstration  has  been  dependent  upon 
the  Initial  conditions  used  from  the  VSLNQH  code  and  the  Kang  and  Dunn 
homogeneous  chemical  kinetics.  Even  though  improved  flow  characterization 
and  verification  are  required,  this  two-dimensional  boundary-layer  code 
has  demonstrated  how  the  oxygen  recombination  on  silicon  dioxide  model 
can  be  used. 


Figure  18.  Mass  Fraction  Profiles  at  x 
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VII.  Conclusions  and  Recommendations 


This  research  dealt  with  the  building  of  a  detailed,  steady-state, 
theoretical  model  for  the  recombination  of  atomic  oxygen  on  a  silicon 
dioxide  surface.  The  Langmuir-Rideal  heterogeneous  recombination  mechani 
formed  the  basis  for  the  model.  The  bonding  of  atomic  oxygen  to  the 
surface  and  the  thermal  desorption  of  atomic  oxygen  from  the  surface 
were  also  included  in  the  model.  The  hypothesis  was  made  that  the  gas- 
phase  oxygen  atoms  combine  directly  with  the  oxygen  atoms  that  constitute 
the  silicon  dioxide  surface,  with  other  gas-phase  oxygen  atoms  replacing 
the  lost  atoms  that  were  bonded  on  the  silicon  dioxide  surface.  A  set 
of  two-dimensional,  steady-state,  laminar  boundary-layer  calculations 
was  made  using  numerical  methods  to  demonstrate  the  use  of  this  model 
and  to  explore  the  rational  limits  which  may  be  placed  on  the  role  of 
catalysis  for  oxygen  recombination  on  a  silicon  dioxide  surface. 

The  model  agrees  with  the  experimental  data  that  is  available  in 
the  literature,  and  provides  an  insight  into  the  processes  that  control 
the  recombination  as  a  function  of  temperature.  From  300°K  to  740°K 
where  0,  the  fractional  surface  concentration  of  atomic  oxygen,  is 
approximately  unity  due  to  rapid  adsorption  of  atoms  to  the  surface,  the 
recombination  process  is  first  order  with  respect  to  the  pressure  of  the 
oxygen  atoms  and  is  controlled  by  the  recombination  desorption  rate.  At 
more  elevated  temperatures  (740°K  to  approximately  1600°K)  where  0  is 
smaller,  the  recombination  process  is  controlled  by  the  chemisorption 
rate  and  remains  first  order.  At  high  temperatures  (above  1600°K)  where 
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the  thermal  desorption  of  atomic  oxygen  is  high,  0  Is  small,  and  the 
recombination  process  Is  second  order  and  is  dependent  upon  the  chemisorp¬ 
tion  rate,  the  thermal  desorption  rate,  and  the  recombination  desorption 
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rate,  A  maximum  value  for  the  recombination  coefficient  Is  1.4  x  10 

at  900°K  which  may  have  a  significance  in  some  applications. 

A  significant  effect  was  found  when  realistic  recombination  rates 

predicted  by  the  model  were  used  in  the  boundary-layer  code.  Although 

the  velocity  profiles  and  temperature  profiles  where  not  appreciably 

changed,  the  mass  fraction  profiles  show  significant  differences. 

In  terms  of  heat  transfer,  the  model  demonstrated  the  most  significant 

Influence.  The  total  heat  load  on  a  reentry  vehicle  Is  highly  dependent 

upon  the  catalyclty  of  the  surface.  In  comparison  to  conditions  where 

the  surface  is  fully  catalytic  and  the  total  heating  rate  Is  0.003494 
2 

kcal/cm  -sec  (x  *  89.  cm),  a  noncatalytlc  surface  would  experience  a 

2 

heating  rate  of  0.002626  kcal/cm  -sec,  with  a  silicon  dioxide  surface 

2 

experiencing  0.003177  kcal/cm  -sec,  for  the  case  where  radiation  equilibrium 
is  assumed.  At  these  conditions,  though  the  silicon  dioxide  recombination 
rate  predicted  by  the  model  Is  low.  It  would  be  a  mistake  to  consider  it 
noncatalytlc;  what  can  be  said,  however.  Is  that  the  predicted  catalycity 
of  the  silicon  dioxide  surface  did  reduce  the  overall  heating  by  9.1% 
below  that  which  would  be  predicted  by  a  fully-catalytic  surface. 

It  has  been  established  that  the  Langmuir  Rideal  mechanism  is  capable 
of  describing  heterogeneous  oxygen  recombination  on  silicon  dioxide. 

The  model  Is  a  realistic  method,  based  on  independent  experimental  data, 
for  predicting  the  recombination  rate  which  Is  required  In  flow  programs 
for  analyzing  the  reentry  of  vehicles  into  the  earth's  atmosphere. 


Further  experimental  work  Is  necessary  to  verify  the  very  high  temperature 
recombination  coefficient  results.  Also,  analytical  and  experimental 
work  is  required  in  the  area  of  energy  transfer  to  the  surface  as  a 
result  of  the  recombination  on  the  surface,  and  how  to  properly  treat 
the  species  in  the  gas  stream  that  might  leave  the  surface  in  an  excited 
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Hydrogen  Recombination  on  Silicon  Dioxide  Model 

The  heterogeneous  recombination  of  hydrogen  atoms  on  silicon  dioxide 
surfaces  has  been  extensively  studied.  It  Is  well  established  that 
hydrogen  atoms  do  chemisorb  to  silicon  dioxide  and  that  hydrogen  molecules 
do  not  (Ref  55-56).  More  recently,  Flnlayson-Pitts  (Ref  57)  confirmed 
that  hydrogen  atoms  chemisorb  to  silica  forming  stable  bonds  and  further¬ 
more  that  hydrogen  molecules  do  not  chemisorb  to  silica— neither  as 
molecules  nor  dissoclatlvely. 

Figure  20  shows  the  recombination  coefficient  as  a  function  of 
temperature  from  three  Independent  sources.  The  data  points  shown  by 
the  stars  of  David  are  from  Smith  (Ref  58)  as  corrected  by  Wood  and  Wise 
(Ref  19).  Wood  and  Wise  determined  that  the  order  of  the  recombination 
process  changes  from  first  to  second  at  a  temperature  greater  than  500°K, 
and  their  data  are  shown  as  triangles  on  Figure  20.  The  data  shown  at 
temperatures  greater  than  500° K  are  the  maximum  values  for  the  recombina¬ 
tion  coefficient  due  to  their  assumption  that  all  of  the  energy  liberated 
upon  the  recombination  of  atoms  went  Into  the  surface.  The  third  set  of 
data  Is  shown  as  squares  and  is  In  the  temperature  range  where  the  recombina 
tion  mechanism  Is  first  order  (Ref  59).  In  contrast  to  the  data  for 
oxygen  recombination  on  silicon  dioxide,  the  data  for  hydrogen  recombination 
from  the  three  different  groups  of  experimentalists  differ  by  less  than 
a  factor  of  three.  Unfortunately,  only  Wood  and  Wise  performed  their 
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Figure  20.  Hydrogen  Recombination  Coefficient  Data  vs.  1/Temperature 
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experiments  at  elevated  enough  temperatures  to  find  where  the  recombination 
kinetics  changed  from  first  to  second  order. 

Wood  and  Wise  (Ref  19)  proposed  that  a  Langmulr-Rldeal  mechanism 
dominates  In  the  temperature  range  300°K  to  500°K,  and  that  the  recombina¬ 
tion  mechanism  Is  the  Langmuir-Hinchelwood  (recombination  of  two  adjacent, 
chemisorbed  atoms)  at  the  elevated  temperatures  where  the  recombination 
kinetics  are  second  order.  However,  the  Langmulr-Rldeal  recombination 
mechanism  can  be  used  to  model  the  recombination  process  for  hydrogen 
atoms  on  silicon  dioxide  above  300°K. 


Analytical  Form  of  Hydrogen  Recombination  Coefficient 

The  exact  same  components  of  the  rate  equation  that  were  described 
In  Chapter  III  are  appropriate  for  hydrogen  recombination  as  well.  For 
steady  state  conditions,  the  recombination  coefficient  can  be  given  by: 


2  P  S0  N  exp(-  T-y) 


B' 


SQ  N  +  5  +  P  N  exp(-  -py) 

B 


(52) 


For  0,  the  fractional  surface  concentration,  to  be  approximately  unity 
at  room  temperature,  the  SQ  N  term  In  the  denominator  must  dominate. 
The  recombination  coefficient  Is  then  given  by: 


Y  =  2  P  exp(-  "j—y)  (53) 

B 

and  the  recombination  kinetics  are  first  order.  Since  the  empirical 
data  from  Wood  and  Wise  show  the  kinetics  going  second  order  above  500°K, 


6,  the  thermal  desorption  rate  per  unit  area,  must  dominate  the  denominator 


of  equation  (52)  at  these  elevated  temperatures.  In  order  for  this  to 
be  true,  the  energy  required  to  thermally  desorb  an  atom  from  the  silicon 
dioxide  surface  must  be  significantly  less  than  It  is  for  oxygen.  The 

•  r 

P  N  exp(-  -py)  term  will  never  dominate  the  denominator  of  equation  (52). 
B 


Hydrogen  Recombination  Data 

An  Arrhenius  plot  of  the  recombination  coefficient  versus  the 

reciprocal  of  the  absolute  temperature  is  linear  for  hydrogen 

recombination  on  silicon  dioxide  from  room  temperature  to  approximately 

500°K.  From  this  plot,  the  activation  energy  Is  2.10  kcal/mole  (1.46  x 
-20 

10  J/atom)  and  the  steric  factor  is  0.06. 

Johnson  (Ref  60)  and  Langmuir  (Ref  61)  found  that  the  number  of 

hydrogen  atoms  adsorbed  per  square  centimeter  of  silicon  dioxide  surface 

15 

is  approximately  1.35  x  10  . 

The  bond  energy,  0,  of  the  chemisorbed  hydrogen  atoms  onto  silicon 
dioxide  Is  reported  by  Wood  and  Wise  (Ref  19)  to  be  approximately  44 
kcal/mole.  They  state  that  Hlrschfelder's  empirical  rule  (Ref  62)  was 
used  to  make  this  estimate  from  the  experimentally  determined  activation 
energy,  E,  of  the  recombination  process.  Hlrschfelder's  rule  Is  the 
following: 

E  =  0.055  D  (93) 

If  the  above  value  of  2.1  kcal/mole  Is  used  for  the  activation  energy, 
the  chemisorption  bond  strength  would  be  38  kcal/mole  from  this  rule. 

Hlckmott  (Ref  32)  performed  thermal  desorption  studies  of  both 
physically  adsorbed  hydrogen  atoms  and  chemisorbed  hydrogen  atoms  In 
addition  to  hydrogen  recombination  tests.  His  results  can  be  interpreted 


such  that  his  molecular  flow  due  to  chemisorbed  atom  recombination  at 
300° K  Is  due  to  thermal  desorption  of  chemisorbed  hydrogen  atoms  followed 
by  Langmulr-Rldeal  recombination  of  these  thermally  desorbed  atoms  with 
chemisorbed  atoms.  With  this  Interpretation,  the  chemisorption  bond 
strength  would  be  approximately  25  kcal/mole. 

Johnson  (Ref  60)  further  found  and  Shuler  and  Laldler  (Ref  63)  and 
Laldler  (Ref  16)  explain  that  thermal  desorption  of  atomic  hydrogen 
becomes  important  at  about  250°C. 

In  addition  to  measuring  the  recombination  coefficient.  Wood  and 
Wise  (Ref  19)  report  that  the  recombination  process  transitions  from 
first  to  second  order  kinetics  with  respect  to  the  pressure  of  the 
hydrogen  atoms  at  a  temperature  of  approximately  500°K. 

Details  of  Hydrogen  Recombination  Model 

The  values  for  the  parameters  In  the  hydrogen  recombination  model 
are  shown  in  Table  VI.  The  Initial  sticking  coefficient  given  in  Table 
VI  is  not  based  Independently  on  experimental  data.  Antonini  (Ref  37) 
found  that  the  initial  sticking  coefficient  for  oxygen  molecules  and  for 
carbon  dioxide  molecules  chemisorbing  to  fresh  silicon  dioxide  decreases 
with  increasing  temperature.  The  same  exponential  form  Is  used  for 
oxygen  atoms  on  silicon  dioxide,  but  the  pre-exponential  factor  was 
Increased  to  get  a  better  match  with  the  recombination  coefficient  data 
shown  In  Figure  20.  The  final  form  Is  as  follows: 

SQ  =  0.75  exp  (-0.002  T)  (94) 

with  the  initial  sticking  coefficient  decreasing  from  0.4  at  300°K  to 


Table  VI.  Computational  Parameters  for  Hydrogen 
Recombination  on  Silicon  Dioxide 


PARAMETER 
STERIC  FACTOR 
E,  ACTIVATION  ENERGY 

S0,  INITIAL  STICKING  COEFFICIENT 
C  .  SURFACE  SITES 

a 

D,  THERMAL  DESORPTION  ENERGY 


VALUE 

0.06 

2.10  kcal/mole  * 

1.46  x  10"^  J/atom 
0.75  exp (-0.002  T) 

1.35  x  10*5  sites/cm^ 

25  kcal/mole  * 

1.737  x  10'19  J/atom 


Figure  21  gives  the  effect  of  changing  the  Initial  sticking  coefficient 

IQ 

with  the  desorption  energy  fixed  at  25  kcal/mole  (1.737  x  10  J/atom). 

The  primary  influence  of  the  Initial  sticking  coefficient  is  at  the 

elevated  temperatures.  We  observe  from  equation  (52)  that  at  the  low 

temperatures  where  the  SQ  N  term  dominates  the  denominator,  the  initial 

sticking  coefficient  cancels  out.  As  long  as  S0  Is  much  greater  than 

P  exp  (-  -py)  at  the  low  temperatures,  then  gamma  will  be  Independent  of 
B 

the  Initial  sticking  coefficient.  At  an  elevated  temperature,  an  increase 
In  SQ  will  Increase  gamma.  Furthermore,  increasing  the  Initial  sticking 
coefficient  moves  the  maximum  value  of  ganina  to  higher  temperatures. 

Figure  22  shows  the  Influence  of  the  hydrogen  atom  desorption  energy 
on  the  value  of  the  recombination  coefficient  with  the  other  variables 
given  by  Table  VI.  The  recombination  coefficient  is  independent  of  the 
desorption  energy  at  the  low  temperatures.  However,  at  the  elevated 
temperatures  the  desorption  energy  Is  very  important.  Not  only  Is  the 


magnitude  of  the  recombination  coefficient  affected  but  the  thermal 
desorption  term  determines  at  which  temperature  the  recombination  kinetics 
transition  from  first  order  to  second  order.  In  order  to  satisfactorily 
match  the  data  shown  in  Figure  20,  25  kcal/mole  (1.737  x  10’19  J/atom) 

Is  used  with  the  rest  of  the  variables  given  in  Table  VI. 

Figure  23  shows  not  only  the  empirical  data  given  by  Figure  20  but 
also  shows  the  results  of  the  model.  The  model's  results  match  the 
empirical  data  very  well.  The  maximum  value  of  the  recombination  coef- 
ficient  is  1.3  x  10  at  500°K.  The  recombination  kinetics  transition 
from  first  to  second  order  with  respect  to  the  atomic  hydrogen  pressure 
at  approximately  550°K.  At  1120°K  gamma  Is  1.1  x  10"7.  This  value  is 
in  agreement  with  the  value  of  approximately  zero  obtained  by  Wood  and 
Wise  at  1123°K. 

Gelb  and  Kim  Model  for  Hydrogen  Recombination  on  Silicon  Dioxide 

Gelb  and  Kim  (Ref  64)  modeled  hydrogen  atom  recombination  on  silicon 

dioxide  using  the  Langmulr-Rideal  mechanism.  They  used  a  value  of  42 

kcal/mole  for  the  thermal  desorption  energy  in  accordance  with  the 

Hirschfelder  relationship.  As  a  consequence,  they  report  a  maximum 

value  for  the  recombination  coefficient  at  833°K.  From  room  temperature 

to  833°K  the  recombination  kinetics  would  be  first  order,  according  to 

their  model,  with  the  recombination  mechanism  changing  to  second  order 

above  833°K.  Better  agreement  with  empirical  data  is  obtained  if  a 

-19 

value  of  25  kcal/mole  (1.737  x  10’  J/atom)  is  used  for  the  thermal 
desorption  energy.  Figure  24  shows  the  results  of  both  their  model  and 
the  model  from  this  research;  each  with  25  kcal/mole  used  for  the 
desorption  energy. 


GAMMA 


Figure  21.  Hydrogen  Recombination  Coefficient  vs.  1/Temperature 
for  SQ  =  0.15  L  exp(-0.002  T)  for  L  =  1,  3,  5. 
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Figure  23.  Hydrogen  Recombination  Coefficient  vs.  1 
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